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SECTION  1 
INTRODUCTION 

Blunt  probes  [Hale  and  Hoult  (1965);  Hale  (1967);  Hale,  Hoult 
and  Baker  (1968)]  and  Cerdien  condensers  [Pedersen  (1964);  Rose  and 
k'fddel  (1972);  Conley  (1974);  Croskey,  Hale  and  Leiden  (1977); 

Mitchell,  Sagar  and  Olsen  (1977)]  are  presently  being  flown  on  rocket 
and  balloon  systems  to  measure  electrical  conductivity  in  the  middle 
atmosphere.  The  Cerdien  condenser  has  the  additional  capability  of 
being  able  to  measure  ion  mobility  and  charge  number  density  [Pedersen 
(1964);  Croskey  (1976);  Sagar  (1976)].  Both  of  these  experiments  are 
being  utilized  to  study  ionization  process  in  the  middle  atmosphere. 

In  particular,  such  phenomena  as  midlatitude  sunrise  condition  [Mitchell 
et.  al_ • (1977)]  and  the  D-region  "v/inter  anomaly"  [Mitchell,  Hale,  Olsen, 
Randhawa  and  Rubio  (1972);  Mitchell  and  Hale  (1973)],  a solar  eclipse 
[Baker  and  Hale  (1970)],  a polar  cap  absorption  event  [Hale  (1974)]  and 
the  high-latitude,  middle  atmosphere  during  geomagnetically  disturbed 
conditions  [Olsen,  Mitchell  and  Croskey  (1976)]  have  been  studied  using 
rocket  instruments.  A balloon-borne  blunt  probe  experiment  on  the  recent 
STRATCOM  flights  has  also  proven  useful  in  studying  the  temperature  de- 
pendence and  altitude  dependence  of  electrical  conductivity  in  the  strat- 
osphere [Mitchell  and  Hale  (1973);  Mitchell,  Hale  and  Croskey  (1977)]. 

kith  the  recent  miniaturization  of  the  blunt  probe  to  make  it  com- 
patible with  the  super  Loki  meteorological  rocket  system,  and  with  this 
particular  instrument  now  commercially  available  [Olsen  (1977)], 


the  number  of  blunt  probe  rocket  flights  has  increased  appreciably. 
This,  in  turn,  has  resulted  in  the  need  for  computerized  data  pro- 
cessing and  reduction  techniques  to  improve  both  the  speed  and 
accuracy  of  these  tasks. 
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SECTION  2 

THF.  BLUNT  PROBE  AND  GERDIEN  CONDENSER  EXPERIMENTS 
2. 1 rxper  i.'.cnt  Peso:  l pt  ton 

The  subsonic  blunt  probe  [Hale  et  a^.  (1965);  Hale  (1967)]  and 
Gerdien  condenser  [Pedersen  (1964);  Croskcv  (1976);  Sagar  (1976)] 
experiments  measure  electrical  conductivity  and  in  addition,  the 
Cerdlen  condenser  measures  ion  mobility  and  charge  number  density. 
When  launched  using  a rocket,  the  payload  separates  from  the  motor 
at  apogee  (nominally  at  70  to  80  km)  and  descends  to  the  ground  on  a 
parachute.  The  data  are  telemetered  in  flight  back  to  the  receiving 
station  where  the  information  is  recorded  on  magnetic  tape. 

The  blunt  probe  experiment  was  initially  designed  in  the  nid 
19  u0 ' s [Hale  c_t  a l . (1966);  Hale  (1967)].  The  theory  of  charged 
particle  collection  for  this  instrument  indicates  that  to  first  order 
the  collection  current  is  neither  dependent  on  the  descent  velocity 
nor  the  pendulum  motion  of  the  payload  as  it  descends  on  a parachute 
[Hale  et  .U  . (1965);  lloult  ( 1 9 o 5)  ] . 

With  the  recent  development  of  more  stabilised  parachute  systems 
the  possibility  of  a subsonic  Gcrdien  condenser  experiment  for  mea- 
suring ion  mobility  and  charge  number  density  (both  of  which  are  flow 
dependent  parameters)  is  now  feasible.  Thus,  subsonic.  Gerdien  con- 
denser experiments  fer  flying  on  such  stabilized  parachute  systems 
have  been  recently  developed  [Karrokh  (1975);  Croskey  (1976);  Sagar 
(1)76)]  to  measure  those  electrical  parameters. 


The  particular  instruments  for  which  the  computer  data  reduction 
scheme  in  this  report  was  developed  are  flown  cn  standard  meteorolog- 
ical rocket  systems  such  as  the  Areas  ar.d  super  Loki  rockets.  The 
10CQ  MHz  transmitter  and  modulation  system  is  compatible  with  those 
used  by  the  Meteorological  Rocket  Network  (MRN) , thus  making  It 
possible  to  launch  the  instruments  at  any  MRN  rocket  range  I The 
Meteorological  Rocket  Network  Document  11-64  (1065)]. 

2 . 2 Current-Voltage  Ra lot ionahi p s 

The  blunt  probe  uses  a circular  planar  collector  geometry  fer 
charged  particle  collection  (see  Figure  (2-1)).  The  current  of  col- 
lected charged  particles  is  described  by  the  equation  [l'.ale  (1967); 
Mitchell  (1973)]: 


o 


where  r and  P.  are  the  radii  of  the  collector  and  the  outside  of  the 
guard  ring,  respectively,  and  is  either  the  positive  or  negative 
electrical  conductivity  which  is  defined  as  follows: 


X V e "i+ 

i 

> N . eg,  + N e u 
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(2-2) 


(2-3) 


In  the  above  expressions,  N ) represents  the  concentration  of 

positive  (negative)  ions  of  the  ith  species  and  (y^_)  *-ts 
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respective  mobility  value.  The  concentration  and  mobility  of  free 
electrons  are  and  p^,  respectively,  and  e represents  the  magnitude 
of  the  charge  for  an  electron. 

If  the  collector  voltage  is  a ramp  function,  It  is  preferable  to 
determine  electrical  conductivity  by  measuring  (dl  /dV)  and  using  t tie 
following  equation: 


o 


t 


dl 

♦ 

dV 


(2-4) 


This  approach  removes  the  problem  ol  having  tv'  actually  measure  the 
probe's  current  and  voltage1  in  determining  o ( . 

The  collected  charge  particle  current  is  measured  by  an  elec- 
trometer [Zimmerman  (l‘)71)]  housed  inside  the  probe,  ami  the  electro- 
meter's analog  output  signal  is  converted  t o a negative  putse  wave- 
form having  a frequency  proport  tonal  to  the  elect romet or 's  output 
voltage.  The  pulse  waveform  in  turn  modulates  the  transmitter  output. 
Thus,  the  expression  for  o in  Kq.  (2-4)  can  now  be  written  in  the1 
form 


o 


t 


i 1 /iI,\>ata 


i.'-S) 


where  f ( is  the  modulation  frequency  correspomling  tv>  the  in-flight 
measurement  ol  collected  charge  particle  current  l,.  Prior  tv' 
launch,  a high-valued,  precision  resistor  (K^ ) is  connected  between 
the  collecting  and  return  electrodes,  thus  feeding  a calibration  cur- 
rent to  the  electrometer.  The  resulting  telemetered  waveform  t^,^. 
is  received  through  the  telemetry  system  prlot  tv'  Might  and  recorded 


< 


7 


on  magnetic  tape.  In  Eq.  (2-5),  the  value  of  this  calibration 
resistor  Is  anc*  the  slope  of  the  modulated  calibration  wave- 


form during  the  sweep  portion  Is  (df/dt) 


CAL 


For  a Gerdien  condenser,  the  collector  geometry  consists  of  a 
cylindrical  collector  and  an  outer,  concentric  cylindrical  return 
electrode  (see  Figure  (2-2)).  A voltage  waveform  is  swept  between 
these  two  electro>!es  and  the  resulting  current  of  charge  particles 
collected  on  the  inner  electrodes  is  measured.  Using  a similar 
electronics  system  and  preflight  calibration  procedure  [Sagar 
( 1 9 7 *5 ) ] os  for  the  blunt  probe,  the  expression  in  the  Gerdien 
ccr.c  ?naer ' s linear  region  of  operation  is 


ln(ro/ri)  1 


(df/dt) 


DATA 


2-  l 


RCAL  <Jf/Jt)CAL 


(2-6) 


In  the  above  expression,  r , r und  £ are  the  inner  radius  of  the 
return  electrode  and  the  radius  and  length  of  the  collecting  elec- 
trode, respect ivclv.  Again,  R„,t  is  the  resistor  value  used  in 
parallel  with  the  condenser  to  generate  the  preflight  calibration 

current  resulting  in  the  modulating  frequency  f_t.  , 

' CAL 

If  the  flow  through  the  aspirator  can  be  determined,  then  the 
Gerdien  condenser's  cylindrical  electrode  geometry  also  affords  the 
opportunity  for  measuring  the  Ion  mobility  and  charge  number  density. 
The  reduction  of  the  Gerdien  condenser's  current-voltage  response  to 
obtain  the  mobility  Information  requires  the  determination  of  the 
voltages  at  which  the  different  ion  mobility  groups  are  entirely 


V 


! 


: 


Figure  (2-2)  Cerdion  Condenser 


collected  out  of  the  air  sample  passing  through  the  aspirator.  The 
determination  of  these  voltage  values  requires  a further  analysis  of 
the  probe's  1-V  response  and  was  not  considered  In  the  scope  of  this 
research. 


I 


2 . 3 Experimental  Procedure 

As  mentioned  previously,  the  rocket  experiments  are  usually  con- 
ducted at  Meteorological  Rocket  Network  (MRN)  launch  sites.  Prior 
to  launch,  the  instrument  is  operated  in  the  preflight  calibra- 
tion mode  as  discussed  earlier.  The  resulting  transmitted  signal 
is  received  and  recorded  on  magnetic  tape.  An  example  of  this  wave- 
form showing  the  calibration  frequency  versus  time  is  given  in  Figure  (2—3), 

The  telemetry  system  used  for  obtaining  the  in-flight  data  is 
the  same  as  that  used  during  the  preflight  calibration.  In  addition, 
timing  information  is  recorded  on  another  channel  of  the  tape  recorder. 

KhiLe  the  probe  is  descending  on  a parachute,  a ground-based  radar 
system  measures  the  position  and  velocity  of  the  instrument  as  a func- 
tion of  time  and  thus  the  electrical  conductivity  data,  which  are  also 
recorded  as  a function  of  time,  can  later  be  determined  as  a function 
of  altitude. 

The  actual  reduction  of  the  data  waveforms  occurs  at  a later 
time  in  the  laboratory.  Representative  data  waveforms  for  both  the 
blunt  probe  and  the  Gerdien  condenser  are  also  shown  in  Figure  (2-3). 

As  discussed  earlier,  the  elect rita!  conductivity  values  are  propor- 
tional to  the  designated  slopes  of  the  modulated  data  waveforms  and 


I 


i 


I 


thus,  the  data  reduction  procedure  involves  determining  the  slopes 


of  these,  particular  waveforms  [Hale  (1967);  Mitchell  (1973);  Sagar 
(1976)]. 
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SECTION  3 

DATA  REDUCTION  MINICOMPUTER  SYSTEM 


3 . 1 System  Functions 

3.1.1  Introduction 

The  overall  minicomputer  system  performs  two  general  functions 
in  reducing  the  electrical  conductivity  data,  namely,  data  acquisition 
and  data  processing.  The  data  acquisition  procedure  involves  the 
transfer  und  storage  of  data  from  the  originally  recorded  magnetic 
tape  to  a DEC  RK11/RK05  disk.  In  transferring  the  data,  which  ini- 
tially are  a series  of  negative  pulses  in  the  frequency  range  of  0 
to  200  pps,  the  time  period  between  the  leading  edges  of  consecutive 
pulses  are  measured  using  the  Real-Time  Clock  of  the  DEC  LPS11  Labora- 
tory Peripheral  System  and  the  digitized  values  are  stored  on  the  disk. 

The  processing  of  the  data  Involves  such  tasks  as  waveform  seg- 
mentation and  display,  the  removal  of  spurious  noise  from  the  data 
waveform,  and  actually  determining  the  slope,  i.e.,  (df^/dt)^  or 
(if/dt)^,^  of  the  waveform  for  a designated  time  interval. 

A discussion  of  the  DEC  i'DP  11/10  minicomputer  and  LPS11  Labora- 
tory Peripheral  System,  which  are  inherent  to  both  of  these  data  re- 
duction functions,  will  be  discussed  in  later  sections  of  this  section. 
Also,  a user’s  manual  of  this  system  Is  presented  in  Appendix  A. 


3.1.2  Data  Acquisition  System 

The  block  diagram  of  the  data  acquisition  system  is  showni  in 
Figure  (3-1).  The  data  inputs  originate  from  the  tape  transport  unit 
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FigOre  (3-1)  Data  Acquisition  System 
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CHP096Q  Instrumentation  Recorder).  The  output  waveforms  from  the  tape 
transport  unit  are  of  three  general  typea:  calibration  data,  in-flight 
data  and  timing  data  (recorded  simultaneously  with  the  in-flight  data). 
The  calibration  and  in-flight  data  are  recorded  in  the  form  of  a nega- 
tive-pulse modulated  signal  which  Is  the  input  to  Schmitt  trigger  2 
of  the  Real-Time  Clock.  The  time  periods  between  two  consecutive 
pulses  are  measured  by  the  Real-Time  Clock  as  digital  integer  values. 
Also,  the  accumulation  of  digital  Integer  values  is  used  to  determine 
the  time  of  the  data  waveforms  as  referenced  to  the  launch  time  for 
the  payload. 

An  optional  method  for  determining  the  relative  times  corre- 
sponding to  the  data  waveforms  involves  using  the  simultaneously  re- 
corded timing  data  on  the  magnetic  tape.  The  timing  information  is 
demodulated  to  a fixed  rate,  offset  binary  signal.  Since  the  signal 
is  not  compatible  with  the  LPS11  Laboratory  Peripheral  System's  dig- 
ital input  port,  it  is  fed  into  an  analog  input  channel  and  sampled 
using  Schmitt  trigger  1 fired  by  an  external  oscillator  running  at 
about  30  kHz . 

All  the  digitized  data  are  stored  on  the  RK11/RK05  disk  of  which 
the  storage  capacity  is  more  than  1.2  million  16-bit  words, 

3.1.3  Data  Processing  System 

The  data  processing  system  is  used  in  an  off-line  or  delay-time 
sense  after  data  are  initially  stored  as  described  in  Section  3.1,2, 

The  block  diagram  of  the  data  processing  system  is  shown  in  Kigure 
(3-2).  The  Central  Processor  Unit  (CPU)  of  the  minicomputer  performs 
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all  the  numerical  functions  and  calculations  such  as  determining  a 
straight  line  fit  to  the  data  waveform  hy  the  least -squares  method, 
setting  a criteria  for  waveform  segmentation,  etc.  A Tektronix  603 
Storage  Scope  is  used  to  display  the  waveforms  for  determining  the 
regions  which  provide  (df/dt)  , (df  /dt)_.„.  or  (df  /dt)_,_,.  The 
I.A36  DEC  writer  II  is  used  to  input  the  signal  to  control  the  proces- 
sing and  to  print  out  a hard  copy  of  the  results. 

/ 

3.2  PDF  11/10  Minicomputer  System 

The  FDD  11/10  minicomputer  system  [RT-11  System  Reference  Manual 
(1975)]  Includes  a Central  Processor  Unit  (CPU),  a core  memory,  a 
large  number  of  peripheral  devices  and  extensive  software.  It  provides 
the  d.ua  storage,  processing  and  printout  functions.  The  system  com- 
ponents and  peripherals  connect  to  and  communicate  with  each  other  on 
a single  high-speed  bus  known  as  the  UNIBUS.  Address,  data  and  control 
Information  are  sent  along  the  56  lines  of  the  bus.  The  form  of 
common  lent  ion  is  the  same  for  every  device  on  the  UNIBUS.  Each  device, 
including  memory  locations,  processor  registers  and  peripheral  registers, 
is  assigned  an  address. 

The  Central  Processor  Unit  Is  connected  to  the  UNIBUS  as  a sub- 
system. It  performs  arithmetic  and  logic  operations,  instruction  de- 
coding, and  data  transfers  directly  between  the  input /output  (1/0) 
devices  and  memory . 

The  core  memory  is  viewed  as  a series  of  locations,  with  a number 
(address!  assigned  to  each  location.  The  FDP  11/10  memory  is  designed 
to  accommodate  both  16-blt  words  and  8-blt  bytes.  A 16-blt  word  used 
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for  byte  addressing  can  address  a maximum  of  32K.  words.  However, 
addresses  from  0 to  777g  and  the  top  4,0%^  word  locations  have  been 
reserved  by  the  system  for  the  interrupt  and  trap  handling,  processor 
stacks,  general  registers,  and  peripheral  devices  registers,  and  there- 
fore, a maximum  of  2SK  of  core  are  left  to  be  programmed.  However, 
only  16K  words  of  core  have  been  implemented  in  this  minicomputer 
system.  The  amount  of  mass  storage  in  the  minicomputer  is  another 
important  consideration  in  this  system.  The  FDP  11/10  system  haB  u 
RK11/RK05  disk  and  two  TAll  cassette  tape  drives  which  can  be  used  for 
immediate  mass  storage.  The  RX11/RK05  disk  has  a maximum  storage 
capacity  of  over  1.2  million  16-bit  words  per  disk  and  a data  transfer 
speed  of  11. I microseconds  (ns)  per  word.  The  RK11/RK05  disk  is  fast 
enough  and  has  enough  storage  capacity  for  about  one  hour  of  flight  time. 
For  example,  if  the  average  frequency  of  the  flight  data  is  100  Hz,  this 
requires  160,000  16  bit-words  of  storage  capacity  if  every  data  point  is 
to  be  stored.  Thus,  if  10,000  data  points  are  transferred  from  memory 
location  to  the  disk,  it  only  takes  0,111  second.  Also,  if  the  direct 
memory  access  (DMA)  operation  is  used,  it  takes  less  than  1 millisecond 
(ms). 

The  cassette  drive  system  is  too  slow  to  be  used  while  the  data 
acquisition  system  is  running  since  there  is  not  enough  time  for  tape 
positioning  and  data  transfer  from  memory  to  cassette  tape.  After  the 
data  have  been  processed  and  stored  on  a disk,  the  cassettes  could  be 
used  as  a cheaper  form  of  storage  for  backup. 

The  communication  between  the  user  and  the  minicomputei  is  also 


an  Important  consideration.  The  I’DP  I1/1Q  system  has  a LA36  DEC 
writer  II  which  is  used  ns  a console  for  inputting  the  control  data 
and  printing  out  the  necessary  Information.  The  1.A36  DEC  writer  II 
is  loaded  with  many  practical  functional  and  operator  features  such 
as  30  character  per  second  throughput  (accompl ished  by  a 60  Hz  catchup 
mode).  Infinitely  variable  vertical  forms  adjustment,  variable  forms 
width,  and  multi-part  forms  capability. 

3.3  I.PSlt  Eahoratory  Peripheral  System 

The  I.PS11  Laboratory  Peripheral  System  includes  a programmable 
Real-Time  Clock  with  two  Schmitt  triggers,  a Display  Control  with  two 
12— bi t D/A  converters,  and  a 12-hit  A/D  converter.  It  is  a high  per- 
formance, modular  and  real-time  subsystem  that  interfaces  with  the 
PDP  11/10  minicomputer  via  the  UN  I ID'S.  The  flexibility  of  the  svstem 
makes  it  well  suited  for  a variety  of  applications  such  as  data  collec- 
tion, monitoring  and  reduction.  A block  diagram  of  this  system  is 
shown  in  l’igure  (3-3). 

3.3.1  I.PSKW  Programmable  Real-Time  Clock 

The  LPSKW  Real-Time  Clock  offers  several  methods  for  accurately 
measuring  and  counting  intervals  or  events.  A block  diagram  is  shown 
in  Figure  (3-4).  The  clock  can  be  used  to  synchronize  the  central 
processor  to  external  events,  count  external  events,  measure  intervals 
of  time  between  events,  and  provide  interrupts  at  programmable  inter- 
vals. It  can  also  be  used  to  start  the  nnalog-to-diglt.il  converter  by 
means  of  the  overflow  from  the  clock  counter  or  by  the  firing  of  a Schmitt 
trigger.  Many  of  these  operations  can  be  performed  concurrently. 


The  clock,  will  generate  one  of  five  crystal-controlled  frequencies: 

1 MHz,  100  kHz,  10  kHz,  1 kHz  or  100  Hz,  and  operate  in  any  one  of  four 
programmable  modes:  single  interval,  repeated  interval,  external  event 
timing  and  event  counting  from  zero  base.  The  Real-Time  Clock  may 
also  use  an  external  (Schmitt  trigger)  input  or  a line  frequency  input 
as  a time  base. 

Two  Sclimltt  triggers  which  are  included  with  the  Real-Time  Clock 
can  start  and  read  the  clock,  start  the  A/D  converter  and  cause  pro- 
gram interrupts. 

3.3.2  LrSVC  Display  Control 

The  LPSVC  Display  Control  is  used  to  display  data  in  the  form  of 
a 4,096^q  x 4,096^  dot  array  on  the  scope.  The  Display  Control  (see 
Figure  (3-5))  consists  of  an  M7019  Scope  Control  Module  and  an  A625 
Digital-to-Analog  converter  Module  which  must  be  used  with  the  M7015 
Bus  Control.  Under  program  control,  a bright  dot  may  be  produced  at 
any  point  in  this  array,  or  a series  of  these  dot3  may  be  programmed 
to  produce  a graphical  output. 

Output  operations  of  the  Display  Control,  which  may  output  to  either 
an  X/Y  recorder  or  a display  unit,  are  accomplished  by  loading  the 
status  register  and  the  X or  Y register.  Through  use  of  status  regis- 
ter bits,  the  Display  Control,  which  operates  the  Tektronix  603  Storage 
Scope  has  the  capability  of  intensifying  the  contents  of  X or  Y registers, 
indicating  when  the  scope  in  ready  for  intensification,  providing  erase, 
write-through,  and  non-stop  control  functions  for  the  storage  scope, 
and  enabling  interrupts. 


3.3.3  LFSAD-12  Analog-to-Digital  Converter 

The  LPSAD-12  is  a 12-bit  successive  approximation  A/D  converter 
which  can  sample  analog  data  at  specified  rates  and  store  the  equivalent 
digital  value  for  subsequent  processing.  The  block  diagram  of  the  A/D 
converter  Is  shown  in  Figure  (3-6).  It  consists  of  four  functional 
modules:  an  M7018  A/D  Control  module,  an  A804  A/D  Converter  module, 
an  A406  Sample-and-Hold  module  and  an  A407  8-Channel  A/D  Multiplexer 


module . 


Figure  (3-6)  A/D  Converter  Block  Diagram 


SECTION  4 


DAT a REDUCTION  PROGRAMMING 

This  section  discusses  the  implementation  of  the  software  pro- 
gramming; for  reducing  the  electrical  conductivity  data.  All  of  the 
programs  were  written  in  Fortran  IV  and  Assembly  Language  under  the 
RT-11  system.  Flow  charts  and  listings  of  the  programs  are  Riven  in 
Appendices  R and  C,  respectively. 

4 . 1 Pulse  Frequency  anil  Timing  Measurements 

In  determining  the  pulse  frequency  of  the  data  waveform  the  Real- 
Time  Clock  is  used  to  measure  the  time  interval  between  the  leading 
edges  of  consecutive  pulses.  The  programming  of  the  Real-Time  Clock 
uses  the  status  register  and  buffer/preset  register,  which  are  des- 
ignated by  lb-bit  words  and  located  at  the  addresses  of  170404 
and  170406  respectively.  When  the  status  register  is  loaded,  It 
enables  the  counter  to  count  at  a designated  rate;  it  controls 
the  rate  of  the  base  frequency  (100  Hz  to  1 MHz);  it  causes  an  in- 
terrupt if  its  flag  is  set;  and  it  counts  the  event  timing  from  zero 
base  by  starting  when  a pulse  comes  from  Schmitt  trigger  2.  Schmitt 
trigger  2 is  used  to  shape  the  data  waveform  by  converting  the  neg- 
ative going  data  pulses  to  a train  of  negative  pulses  with  the  same 
frequency  and  a pulse  width  of  approximately  2 ps.  Each  pulse  from 
Schmitt  trigger  2 sets  the  interrupt  flag  and  causes  the  transfer  of 
the  contents  of  the  huf for/preset  register  to  a temporarily  reserved 


buffer  in  memory,  until  all  the  signals  on  magnetic  tape  have  been 
digitized  and  stored  in  digital  form  on  the  disk.  Therefore,  the 
number  of  counts  is  proportional  to  the  time  interval  of  two  con- 
secutive pulses  and  inversely  proportional  to  the  pulse  frequency. 

For  example,  if  the  100  kHz  clock,  frequency  is  used,  the  number  of 

counts  for  a 100  pulses  per  sec  (pps)  waveform  Is  1,000  and  for  a 
250  pps  waveform,  it  is  400, 

However,  it  the  time  interval  of  two  consecutive  pulses  is 
greater  than  0.65535  sec  and  the  base  frequency  of  the  clock  is  100 
kHz,  the  number  of  counts  will  be  greater  than  65535^  which  is  the 
maximum  value  that  the  15- hit  buffer/preset  register  can  handle. 

An  interrupt  service  routine  and  an  Interrupt  waiting  loop  are  used 
when  a time  interval  between  two  consecutive  pulses  is  greater  than 
0,65535  sec.  An  index  number  within  the  interrupt  waiting  loop  Is  set 
proportional  to  the  number  of  times  that  the  waiting  loop  has  been 
completed  between  the  two  consecutive  pulses.  Also,  the  instruction 
time  ot  the  waiting  loop  can  be  calculated.  Thus,  an  estimate  of  the 
waiting  loop’s  total  execution  time  actually  gives  the  time  interval  of 
the  two  consecutive  pulses.  For  example,  if  the  instruction  time  for 
a complete  waiting  loop  is  23,2il0%  ps  and  the  waiting  loop  has  been 
completed  28,248^  times,  the  Index  number  will  decrease  by  1.  Thus, 

If  the  index  number  is  -100,  the  time  interval  of  these  two  pulses  is 
65.535  seconds. 

The  "WRITE"  request  is  a "programmed  request"  that  is  an  assem- 
bler macro  call  written  into  the  program  and  interpreted  by  the 
FDF  11/10’s  monitor  at  the  program  execution  time.  It  is  used  to 
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transfer  a specific  number  of  data  points  from  a temporarily  reserved 
buffer  of  memory  to  the  disk.  The  control  of  the  Central  Processor 
Unit  returns  to  the  program  immediately  after  the  request  is  queued 
(cl  ms).  The  storage  of  data  points  requires  double-buffered  1/0 
techniques,  i.e.,  the  contents  of  one  of  the  buffers  are  transferred 
from  memory  to  the  disk  while  the  other  buffer  is  filled  immediately 
without  interacting  with  the  previous  buffer.  Since  the  data  are 
digitised  and  collected  from  the  tape  recorder  which  is  continuously 
running,  the  "WRITE”  request  interrupts  the  Central  Processor  Unit 
for  less  than  1 ms,  which  means  that  only  pulses  greater  than  1 kHr 
will  be  lost. 

Since  each  block  of  the  DEC  RK11/RK05  disk  contains  40CL  (256,,.) 

o lu 

16-bit  words,  the  number  of  words  of  every  buffer  reserved  in  memory 
is  usually  an  integer  multiple  of  256^,  i.e.,  every  buffer  could 
contain  14400g  (6400^q)  16-bit  words  (25  blocks) . 

The  "TTY IN"  request  Is  used  to  receive  the  characters  from  the  LA36 

DEC  writer  II.  All  of  the  characters  received  are  in  the  form  of  ASCII 

Code.  Thus,  a subroutine  "TTIN"  is  used  to  convert  the  characters  from 

ASCII  Code  to  numerical  values.  For  example,  the  16-bit  word  34465„ 

o 

(the  higher  byte  is  071^  and  lower  byte  065^)  In  ASCII  Code  converts 
to  a numerical  value  of  59jq. 

Typically,  the  data  recorded  on  magnetic  tape  for  one  hour  occupy 
about  900  blocks  (230,400  16-blt  words)  on  the  disk,  of  which  approxi- 
mately 200  blocks  are  for  the  proflight  calibration  waveforms,  80 
blocks  for  the  data  between  the  launch  of  the  rocket  and  the  separa- 
tion of  the  payload  from  the  rocket,  and  the  remaining  620  blocks  for 
the  in-flight  data. 
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4 . 2 Waveform  Segment at  ton  and  Reatorage 

For  an  easy  access  to  each  waveform  for  data  processing,  the  data 
are  segmented  on  an  Individual  waveform  basis  with  the  data  points 
stored  in  seven  blocks  (per  waveform).  The  restored  data  are  assigned 
a different  file  name  on  the  disk,  and  the  original  data  file  is  used 
as  a back-up. 

For  a typical  experiment,  usually  ten  to  twenty  preflight  call- 

| 

I brat  ion  waveforms  and  one  hundred  or  more  in-flight  data  waveforms 

are  recorded  on  n magnetic  tape.  Referring  to  Figure  (2-3),  most 
waveforms  contain  three  portions:  a beginning  portion  of  low  fre- 
quencies, a region  of  increasing  frequencies,  and  a final  portion  of 
relatively  high  and  almost  constant  frequencies.  Segmenting  the  data 
into  individual  waveforms  was  accomplished  by  recognising  the  general 
features  of  the  beginning  portion  of  the  next  waveform,'  In  doing  this, 
several  hundred  to  a thousand  or  more  data  points  are  extracted  for  a 
complete  calibration  or  in-flight  waveform. 

, 4.3  Pat,*  Display 

Although  each  data  waveform  has  certain  general  characteristics, 
the  large  variability  of  the  central  portion  of  the  waveform  leads  to 
no  simple  and  dependable  method  of  software  programming  to  extract  the 
linear  regions  which  will  provide  the  values  of  (df^/  dt)^^^  (from 
which  the  electrical  conductivity  is  derived).  An  example  of  how  a 
complete  waveform  is  displayed  is  shown  in  Figure  (4-1).  The  values  of 
the  horizontal  coordinate  indicate  the  sequential  number  of  the  data 
points  and  the  scale  Is  a function  of  time.  The  values  and  scale  of 
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the  vertical  coordinate  are  .1  function  of  frequency  in  Hz.  There- 
tore,  the  scope  actually  displays  frequency  as  .1  function  of  time. 

"’•o  ic'pc  —.it  r i X is  a 4,(96^  :t  dot  array.  Frequency 

..nd  t "a  .'.it  i c.  gs  c the  waveform  .ire  required  before  displaying  on 
tin*  seep.  . : e rel  » t i . .'.ship  between  the  original  values  of  the 

dat...  u.  cy,  tie.  , and  or  responding  va’.us  on  the  scope  are 

• h.?vn  i t . »h!e  l,  where  I.t'SVi  a represents  the  values  of  the  horizontal 
coordinate  and.  l.TSW  V represents  the  values  of  the  vertical  coordi- 
nate . 


Sine,  the  date  points  which  provide  for  electrical  conductivity 
i”.*'ovr  .-in  are  within  a United  interval  01  the  overall  wave  fern, 
it  i . often  do-irabl  • to  expand  a specific  region  of  the  waveform 
on  C : c scope  in  order  to  better  understand  and  more  accurately 
choose  the  important  data  points.  For  example,  referring  to  Figure 
\ *-l)  , it  a region  from  the  200rh  to  *00th  data  points  of  the  hori- 
zontal coord  innt  e and  M>  to  ICO  (I.:  ci  the  vertical  coordinate  are 
redispl  lye’.  the  resulting  expanded  waveform  occurs  as  shown  in 
Figure  v«*--2).  The  relationship  of  the  original  values  of  data, 
frequency,  time  and  the  corresponding  values  on  the  scope  are  pio- 
vidod  ir.  r.ihlo  2.  ld.entifie.it!  m ef  .1  waveform  segment  requires 
n umber In ; ill  of  the  data  points  in  sequential  order.  The  numerical 
character.?  (from  0 to  ) are  also  displayed  on  the  scope  and  are 
constructed  using  discrete  dots  on  a S x 7 grid  as  shown  in  Figure 


(s-  n . 


33 


4.4  Data  Processing 

According  to  Eqs.  (.2-5)  and  (.2-6),  the  electrical  conductivity 
information  is  obtained  from  the  values  of  (df/dt)^^  and  (df+/dt) 
DATA  are  c^e  sl°Pes  of  specific  regions  of  the  calibration 

and  data  waveforms,  respectively.  The  slopes  are  obtained  by  fit- 
ting a one-independent  variable  equation,  i.e.,  a straight  line,  to 
the  designated  data  points  of  the  waveform.  The  equation  is  ex- 
pressed as  follows: 

f = A t + B (4-1) 

where  f is  the  frequency  in  Hz,  t is  the  time  in  seconds,  A is  the 
slope  of  the  straight  line  in  Hz/second  (which  is  also  the  value  of 
(df/dt))  and  B is  the  intercept  in  Hz.  The  straight  line  generated 
by  this  equation  is  also  displayed  on  the  scope  to  assure  that  it 
fits  within  the  specified  tolerance.  Special  techniques  for  deter- 
mining a least-squares  straight  line  fit  to  the  data  are  presented 
in  Section  5, 

4.5  Timing  Data 

The  related  real  time  (referenced  to  the  payload's  launch)  at 
which  a data  waveform  is  scaled  can  be  obtained  by  accumulating  the 
interval  for  all  of  the  data  pulses.  This  time  information  is  impor- 
tant in  later  determining  the  payload’s  altitude  from  radar  data, 
which  are  also  usually  referenced  to  launch  time. 

There  are  three  possible  cases  to  consider  in  converting  the 
number  of  counts  of  data  to  the  value  of  real  time.  If  the  number 
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of  counts  of  a data  point  is  less  than  32,767^^  (077777^),  the  real 
time  CTr)  fur  this  data  point  Is  simply  accumulated  as  follows; 


T 

r 


N 

100,000 


+ T. 


(4-2) 


where  is  the  real  time  accumulated  before  this  data  point,  N is 
o 

the  number  of  the  counts  of  this  data  point.  If  the  number  of  counts 
is  greater  than  32,767^  (077777g)  and  less  than  65,535jq  (177777 g) t 
the  real  time  for  this  data  point  lias  to  be  expressed  a*.  . 


(65,536  + N) 
100,000  b 


(4-3) 


If  the  number  of  counts  is  greater  than  65,535^  (177777g),  the 
data  acquisition  system  generates  a particular  v. lue  (-IOO^q) , which 
is  stored  before  the  index  number  in  order  to  bt  easily  recognized 
and  the  index  number  "I"  is  transferred  to  the  real  time  as  follows; 


Tf  - 0.65535j0  X (-1.  + Tg  (4-4) 

where  the  index  number  I is  stored  in  the  data  acquisition  system 
as  a negative  value.  The  time  determined  by  accumulating  the  data 
counts  is  estimated  to  be  within  1%  of  the  actual  time  value,  which 
is  considered  satisfactory  for  a subsonic  experiment. 


SECTION  5 


CURVE  FITTING  TECHNIQUES 

This  section  discusses  a method  of  fitting  a straight  line 
equation  to  a mass  of  data  [Daniel  and  Wood  (1971)].  The  method  for 
determining  a straight  line  fit  should  use  all  of  the  relevant  data 
in  estimating  each  constant  of  the  equation,  have  reasonable  economy 
in  the  number  of  constants  required,  provide  6ome  estimate  of  the  un- 
controlled error  and  give  some  idea  of  how  well  the  final  equation 
can  be  expected  to  predict  the  response. 

5 . 1 Linear  Least-Squares  Method 

The  most  popular  method  of  fitting  an  equation  to  a mass  of  data 
Is  the  least-squares  method.  This  method  finds  the  values  of  the 
constants  in  the  designated  equation  such  that  the  sum  of  the  squared 
deviations  of  the  observed  values  from  those  predicted  by  the  equation 
is  minimized. 

From  the  data  pairs  (tnl,  fnl>,  (tnl+1,  fnl+]) (tn2,  fn2>  of 

the  specified  region  of  the  waveform  (see  Figure  (5-1)),  there  are  four 
assumptions  about  the  relationship  between  the  observed  value  of  the 
Independent  variable  t^  and  the  observed  value  of  the  dependent  vari- 
able f^  for  determining  a straight  line  y = A x + B to  fit  the  data: 

1.  There  is  a linear  relationship  between  the  predicted  value  of 
a response  y and  the  value  of  the  independent  variable  x 

y - A x + B (5-1) 


where  A is  the  slope  and  B is  the  intercept. 
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2.  The  observed  value  ■ A + B + e^,  where  Is 

the  random  error. 

3.  The  random  error  e^  has  the  following  properties: 

1)  The  expected  value  of  e^  is  zero  and  the  observed  f is 
an  unbiased  estimate  of  y . 

2)  The  variance  of  e^  remains  constant  for  all  values  of  x^. 

3)  The  Oj  values  are  statistically  uncorrelatcd. 

4.  The  observed  values  of  the  independent  variable  are  measured 
without  error.  All  of  the  error  is  thus  in  f^  and  none  is  in 
the  t^'s. 

Therefore  the  linear  least-squares  estimates  for  a straight  line 
are  those  values  of  A and  N which  minimize  the  function  as  follows: 


n2  n2 

n2 

q - £ c^)-’  - £ ■ 

X <trAirB)2 

(5-2) 

i“nl  i=*nl 

i^nl 

From  the  above  assumptions,  there  is  no  error  for  the  Independent 
variable  t^  and  hence,  x^  **  t^.  Thus  we  obtain 
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Thus,  the  predicted  value  Is  a function  of  the  known  data  pairs 
(tj,  f^)  only  and  can  be  expressed  as  follows: 

y - A x£  + B - A t + B 
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Finally  the  straight  lino  y “ A x + B is  derived  in  which  the  slope 
A is  the  value  of  (J f/dt >CAL , (df+/dt)DATA  or  (df_/dt )DATA  tn  *q- 

(2-5)  or  (2-6). 

5.2  Confidence  Rep, ion 

The  digital  data  acquired  from  magnetic  tape  in  general  contain 
some  unexpected  noisy  data.  Therefore  a confidence  region  is  defined 
by  an  appropriate  factor  which  corresponds  to  the  error  between  the 
observed  value  f^  and  the  predicted  value  y^  of  the  fitted  equation. 
According  to  Figure  (5-1),  the  confidence  region  is  between  the  two 
straight  lines  LI  and  L2  for  which  every  point  can  be  expressed  as 
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follows : 


Ll: 

L2 : v. 


*i  + eo 
yi  - “o 


(5-11) 

(5-12) 


where  o is  the  value  chosen  to  define  the  confidence  region.  For 
o 

the  jth  data  point  within  the  specified  region  of  the  waveform,  the 
error  e.  between  the  predicted  and  observed  values  is 


o n V — f 

J ’I  J 


(5-13) 


If  Oj  > e^ , obviously  the  jth  data  point  (such  as  (t^»  f^)  in 
Figure  (5-1))  is  classified  as  a noisy  data  point  and  is  not  included 
in  Fqs . (5-2)  through  (5-10).  Therefore,  the  error  caused  by  the 
unexpected  noisy  data  point  is  eliminated. 


5 . 3 Residual  Root  Moan  Square 

The  residual  root  mean  square  is  computed  to  determine  how  well 
the  final  equation  can  be  expected  to  predict  the  response.  The  value 
of  the  residual  root  mean  square  is  expressed  as 
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where  the  jlth,  J2th,...  data  points  are  not  included  in  Eqs.  (5-2) 

through  (5-10)  *o  obtain  A and  n of  Fq.  (5-1).  The  final  number  of 

total  data  noints  N of  Fqs.  (5-R)  and  (5-*>)  is  always  less  than  or  equal 
o 

to  N. 


SECTION  6 


DISCUSSION  AND  CONCLUSIONS 

6.1  Electrical  Conductivity  Measurements 

The  computerized  reduction  system,  as  described  in  the  previous 
sections,  was  used  to  reduce  blunt  probe  electrical  conductivity  data 
obtained  from  a rocket  experiment  launched  at  1230  MST  on  September  28, 
1976  from  White  Sands  Missile  Range,  New  Mexico.  This  particular  rocket 
experiment  was  conducted  in  conjunction  with  the  STRATCOM  VILA  balloon 
flight  launched  from  Holloman  Air  Force  Base,  New  Mexico.  The  dots  in 
Figure  (6-1)  represent  the  conductivity  values  obtained  using  the  re- 
duction program.  Smooth  curves  have  been  fitted  to  the  data. 

Abcve  40  km,  the  profiles  diverge  with  the  negative  conductivity 
values  larger  than  the  corresponding  positive  conductivity  measurements. 
In  this  region,  the  negative  charge  particles  are  generally  more  mobile 
than  the  positive  Ions.  Below  40  km,  the  positive  and  negative  conducti- 
vity values  for  the  same  altitude  are  comparable. 

The  pulses  in  Figure  (6-1)  represent  conductivity  values  resulting 
from  manually  scaling  the  demodulated  waveforms  (see  Figure  (2-3))  to 

obtain  the  values  of  (df  /dt)_._.  and  (df/dt)_,, . This  is  the  procedure 

i DATA  CAL 

formerly  used  for  reducing  electrical  conductivity  data  [Hale  and  Hoult 
(1965);  Mitchell  (1973)]. 

It  Is  also  important  to  note  that  the  computer  reduction  method 
has  expanded  the  altitude  region  over  which  the  conduct ivity  data  wave- 
forms can  be  reduced.  At  higher  altitudes  where  the  slopes  of  the  data 
waveforms  are  the  largest,  it  is  very  difficult  to  manually  scale  them 


Figure  (6-1)  Parachute-Borne  Blunt  Probe  Electrical  Conductivity  Data 
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from  the  strip  chart.  Also,  the  response  times  of  the  data  demodulation 
system  and  the  strip  chart  recorder  can  introduce  possible  errors.  At 
lover  altitudes  where  the  data  waveforms'  slopes  are  very  small,  It  is 
also  very  difficult  to  manually  scale  the  data.  By  using  the  computer- 
ized scheme  to  reduce  the  data  for  the  September  28,  1976,  blunt  probe 
experiment,  values  were  obtained  for  the  electrical  conductivity  In 
the  25  to  70  km  altitude  region  where  this  parameter  is  observed  to 
change  by  approximately  four  orders  of  magnitude. 

Another  check  of  this  reduction  scheme  is  to  compare  the  results 
of  the  rocket  flight  with  the  data  from  the  STRATCOM  VIIA  blunt  probe 
experiment,  which  are  shown  in  Figure  (6-2).  Included  in  this  figure 
are  the  smooth  curves  for  the  rocket  data  in  Figure  (6-1)  and  the  bal- 
loon blurt  probe  values  obtained  from  1100  to  1800  MST  on  September 
28,  1976,  while  the  balloon  slowly  descended  from  39  km  to  19  km.  The 
balloon  data  format  made  it  more  suitable  for  reduction  by  manually 
scaling  the  waveforms  fron  a strip  chart.  Again,  very  good  agreement 
was  observed  between  the  two  sets  of  data.  This  hopefully  suggests 
consistency  between  both  the  data  reduction  procedures  and  the  exper- 
iment techniques. 

6 . 2 Cone  1 us  tor s 

A computerized  system  using  the  DEO  PDF  11/10  minicomputer  and 
associated  peripherals  has  been  developed  for  reducing  subsonic  blunt 
probe  and  Oerdien  condenser  electrical  conductivity  data.  Assembly 
Language  and  FORTRAN  IV  programs  were  written  under  the  PEC  RT-11 
operating  svstrn  to  perform,  data  digitizing,  acquisition,  storage, 
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display,  processing  and  printing  out  the  results.  Interaction  with  the 
system  is  necessary  for  interpreting  the  data  waveforms  and  for  choosir. 
the  particular  segments  of  the  waveforms  to  be  scaled. 

The  electrical  conductivity  values  determined  by  using  this  tech- 
nique on  rocket  blunt  probe  data  were  found  to  be  consistent  with  the 
Values  obtained  by  manually  scaling  the  data  waveforms.  In  addition, 
the  enhanced  accuracy  of  the  computerized  system  makes  it  possible  to 
scale  the  data  over  a larger  altitude  region,  i.e.,  a wider  range  of 
values,  than  was  found  possible  using  the  manual  technique. 

In  the  future,  a noninteractive  system  using  a microprocessor 
with  a minicomputer  on-line  to  handle  the  data  acquisition  and  process- 
ing simultaneously  would  provide  a substantial  savings  in  the  time 
required  for  data  reduction.  Implementation  of  such  a system,  however, 
would  require  considerable  attention  to  developing  the  computer  soft- 
ware for  interpreting  the  probe's  current-voltage  response. 
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APPENDIX  A 

USER’S  MANUAL 

This  section  was  written  to  assist  the  user  In  Interfacing  the 
system,  using  the  software  programs  and  obtaining  the  results. 

Al.  Data  Acquisition  System 

1)  Equipment 

1.1  DEC  PDP  11/10  Minicomputer 

1.2  Real-Time  Clock  of  LPS11  Laboratory  Peripheral  System 

1.3  HP3960  Instrumentation  Recorder  (or  equivalent) 

1.4  LA36  DEC  writer  II 

1.5  Dual  trace  Oscilloscope 

2)  Instructions 

2.1  Connect  the  output  terminal  of  the  data  channel  of  the  HP396C 
Instrumentation  Recorder  to  both  of  the  input  terminals  of 
Schmitt  trigger  2 of  the  LPS11  Real-Time  Clock  and  channel 

1 of  the  dual  trace  oscilloscope. 

2.2  Connect  the  output  terminal  of  Schmitt  trigger  2 of  the 
LPS11  Real-Time  Clock  to  the  input  terminal  of  channel  2 of 
the  dual  trace  oscilloscope. 

2.3  Set  "-Slope"  of  Schmitt  trigger  2. 

2.4  Adjust  the  level  of  Schmitt  trigger  2 until  it  fires  in  the 
correct  position  ns  observed  on  the  oscilloscope, 

3)  Software  Program  Usage 

3.1  Input  "R  NEW102"  to  the  console  of  the  LA36  DEC  writer  II. 

3.2  Input  the  value  of  the  number  of  blocks  for  every  buffer  to 


the  console. 
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3.3  Input  the  values  of  the  times  of  transferring  data  from 
memory  to  the  disk  to  the  console. 

A2.  Waveform  Segmentation  and  Restorage 

1)  Equipment 

1.1  DEC  PDP  11/10  Minicomputer 

1.2  Display  Control  of  LPS11  Laboratory  Peripheral  Svstem 

1.3  Tektronix  603  Storage  Scope 

1.4  LA36  DEC  writer  II 

2)  Software  Program  Usage 

2.1  Input  "R  NEW500"  to  the  console  of  the  LA36  DEC  writer  II. 

2.2  Find  the  data  at  the  launching  of  the  rocket.  Input  "3" 
to  the  console. 

2.3  Find  the  data  between  the  launching  and  the  separation  of  the 
payload  from  the  rocket.  Input  "4"  to  the  console. 

2.4  Find  the  data  at  the  separation  of  the  payload  from  the 
rocket.  Input  "5"  to  the  console. 

2.5  Input  the  values  of  the  number  of  blocks  for  every  buffer, 
the  total  number  of  blocks  of  the  input  data  file  and  the 
total  number  of  segmented  waveforms,  respectively,  to  the 
console . 

A3.  Data  Processing 

1)  Equipment 

1.1  DEC  PDP  11/10  Minicomputer 

1.2  Display  Control  of  LTSll  Laboratory  Peripheral  System 

1.3  Tektronix  603  Storage  Scope 

1.4  LA36  DEC  writer  II 
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2}  Software  Program  Usage 

2.1  Input  "R  NEW900"  to  the  console  of  the  LA36  DEC  writer  II. 

2.2  Input  the  sequential  number  of  the  first  waveform  to  be 
processed  to  the  console. 

2.3  Input  "1"  to  the  console  for  expansion  of  the  specific  region 
of  the  waveform. 

2.4  Input  "2"  to  the  console  for  obtaining  the  positive  electrical 
conductivity. 

2.5  Input  "3"  to  the  console  for  obtaining  the  negative  electrical 
conductivity. 

2.6  Input  "4"  to  the  console  for  the  processing  of  the  other 
waveform. 

2.7  Input  ''5"  for  terminating  the  data  processing  and  printing  out 
the  results. 

2.8  Input  the  sequential  numbers  of  the  first  and  last  data  points 
and  the  tolerance  of  the  frequency  to  the  console  for  obtaining 
the  electrical  conductivity. 
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APPENDIX  B 


FLOW  CHARTS  OF  MINICOMPUTER  (RT-11)  PROCRAMS 
B1  - Flow  chart  of  data  acquisition. 

B2  - Flow  chart  of  waveform  segmentation  and  restorage. 

B3  - Flow  chart  of  data  processing  (using  a least-squares  method) 
and  printing  out  the  results  . 
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CALL  CMPR1 


Look  for  the  beginning  part  of  the 
waveform,  reserve  data  points  in 
buffer,  input  data  from  disk  if  nee- 
essarv 


CALL  CMPR 


CALL  CMPKJ 


Write  data 
into  Disk 


Fm  i sli? 


Look  for  the  central  part  of  the  wave- 
form, reserve  data  points  in  buffer, 
input  data  from  disk  if  necessary 


Look  for  the  beginning  part  of  the 
next  waveform,  reserve  data  points 
buffer,  input  data  from  disk  if  no 


Flow  Chart  of  Waveform  Segmentation  and  Restorage 
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CALL  SL0P1 


L'se  least-square  method  to  find 
the  slope  of  the  specified  linear 
region  of  the  waveform  and  the 
designated  time  interval,  remove 
the  spurious  noise,  display  the 
fit  straight  line  on  scope 


CALL  SLOP 


Use 

least-sqi 

are  method  to 

find 

the 

slope  of 

the  specified 

1 inear 

region  of  the 

waveform  and 

the 

des] 

ignated  time  interval,  . 

•?move 

the 
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fit 

straight 

line  on  scope 

B3.  Flow  Chart  of  Data  Processing 
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APPENDIX  C 

MINICOMPUTER  (RT-11)  PROGRAMS 


This  section  gives  the  source  programs  which  are  run  on  the  PDF 
11/10  minicomputer.  These  programs  are  constructed  either  by  link- 
ing groups  of  Assembly  Language  programs  or  by  linking  Fortran 
programs  with  Assembly  Language  programs.  The  object  modules  of  the 
main  programs  and  subroutines  are  linked  as  shown  in  Appendix  Cl,  and 
all  the  texts  of  the  programs  are  listed  in  Appendices  C2  through 
C23.  The  functions  of  the  source  programs  are  listed  as  follows: 
NEV102  - Frogram  for  digital  data  acquisition. 

NEW500  - Program  for  waveform  segmentation  and  restorage. 

NEW900  - Program  for  data  processing  (using  the  least-squares 


method)  and  for  printing  out  of  the  results. 


appendix  c 

. R LINK 

+NEN102<NEW10O/ NEW101 


R LINK 

*NEW500<NEW200,'C 

♦NEWlOi/C 

+NEW31 Q> NEW320. NEW330, NEW340/C 
♦NEN400,  NEW420. NEM440/C 
+NEM240, NEU250 


. R LINK 

♦NEW 90O< NEWS 02/F/C 
♦NEN810/G  i/C 
♦NEW825/0 . 2/C 
♦NEM345/0  2/C 

*NEW310. N E W 320. NEW" 33- NEW340/0  3/C 
♦NEW315,  NEW325-  NEW335- HEW345- NEU355/0. 


Cl  - Source  Programs 
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lift  IN 


MILT  I 


START 


WRIT 


. CSECT 
. TITLE' 
GLOBL 
LPSCkS* 
LPSF'B2  = 
. MCALL 
MCALL 
. V2.  . 

. REGDEF 
. PRINT 
J SR 
MOV 
MOV 
MOV 
CIR- 
CLE 
CLR 
CLR 
CLR 
MOV 
MOV 
ADD 
INC 
CMP 
PNE 
MOV 

. FETCH 
ENTER 
PRINT 
SUE 
MOV 
MOV 


N F W 100.  MAC 
MAIN,  TT  I N 
170404 
1 7040b 

. . V2  . , . REGDEF, 
PRINT, . EXIT 


F'TCH, . ENTER, . WRITE,.  CLOSE 


#riSG  i 
FC,  TT!  N 
4DATAIN. F 2 
<R2>*. ELOCkl 
<P2>«-.  CVCLE 
PI 
F 2 
RS 
R5 

EL0CK2 

4SERV.  0*724 

»S40,  c« # : 2 e 

#400-  F-: 

R 1 

ELOCkl, F 1 
MULTI 
R‘ , LREC 
4HNDR,  4N  AME 
(♦AREA,  *1,  4 N A M E » 
(♦MSG  2 

EL0CK1. ELOCkT 
♦♦BUFFI,  F:i 
Ri,  R4 


PRINT  OUT  CONTENTS  OF  MSG1 
RECEIVE  OATH  FROM  DECWRITER 
SET  IIP  R2  WITH  ADDRESS 
TO  ELOCkl 
TO  CVCLE 


DATA 

DATA 


, MOVE 
, MOVE 
, R1  = 0 
, R7  = 0 
, R S = O 
, R 5 * 0 

-INITIALIZE  BLOCKS  WITH  0 
, SET  UP  ADDRESS  AND  PRIGRITV 
, LEVEL  OF  INTERRUPT  ROUTINE 
, R S = 2 5 6 * B L 0 C k 1 
■INCREASE  Rl  EV  1 
, R 1 = E L 0 C k 1 ? 

. NO.  REPEAT  ADDING 
, VES.  R 2 = L R E C =256* ELOCkl 
, DEFINE  FILE 
4-1 

.PRINT  OUT  CONTENTS  OF  MSG 2 
, INITIALIZE  BLOCKS 
SET  UP  Rl  WITH  ADDRESS 
i SET  UP  R4  WITH  ADDRESS 
-SET  UP  EXTERNAL  INTERVAL 
, FROM  ZERO  EASE,  FREOUENCV 
, EEING  100  KHZ,  INTERRUPT 
, ENABLE,  CLOCK  ENABLE 


MOV 

41505.  LPSCkS 

-WAITING  LOOP  FOR  INTERRUPT 

INC 

R2 

; SET  UP  INDEX  OF  REAL-TIME 

CMP 

467110,  R 2 

, EV  THE  INSTRUCTION  TIME 

EGT 

WAIT 

, OF  THE  WAITING  LOOP 

DEC 

R5 

, R5  DECREASE  EV  1,  REAL-TIME 

CLR 

P2 

-PASSED  2S  2*2E.24E*E-6  SEC 

CMP 

40,  re- 

•IF PS=0.  BUFFER  IS  FULL 

EL  T 

start 

. NO,  DO  NEXT 

CLR 

LF'SCkS 

, VES, DISABLE  THE  INTERRUPT 

CMP 

F 4,  4EUFF1 

- DOUBLE-BUFFER  METHOD 

ENE 

ALPHA 

. R4  POINTS  TO  BUFF  1 ° 

MOV 

4EUFF2, Rl 

, VES.R1  POINTS  TO  BUFF2 
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BR 

STORE 

f>  L F A A 

MOV 

#BUFrl,Rl 

N0.R1  POINTS  BUFFI 

? T ORE 

AM* 

BLOCK  1 ■ BLOC  EZ 

SET  UP  BLOCK  FOR  THE  FILE 

MOV 

L F EC, r: 

WRITE 

#ARE A • #1,  R 4 . Rj,  CL 

OCT  2 i TRANSFER  DATA  FROM 

CMP 

R4 . #BUFFi 

, MEMORY  TO  DISK 

BNE 

BETA 

R4  POINTS  TO  BUFFI'’ 

MOV 

# E Ur  F Z.  R4 

VES. R4  POINTS  TO  BUFF2 

BR 

GAMMA 

BETA 

MOV 

#BUFF1,R4 

NO. R 4 POINTS  TO  BUFFI 

GAMMA 

PRINT 

#MSGZ 

PRINT  OUT  CONTENTS  OF  MSGS 

DEC 

CVC  IE 

STORE  ONCE 

BLE 

FINISH 

JMP 

START 

SERV 

CMP 

#0.  R5 

INTERRUPT  SERVICE  ROUTINE 

FEQ 

OMEGA 

MOV 

#-100  , (Rl>  + 

SAVE  THE  COUNT  NUMBER  AND 

MOV 

R 5 • >.  R 1 > * > 

TIMING  INDEX  IF  NECESSARY 

SUB 

#z.  r: 

TWO  DATA  PTS  FOR  TIMING  INDEX 

CLR 

R5 

RESTORE  R5  TO  0 

OMEGA 

MOV 

LFSRBZ. < R1 >♦  ; 

SAVE  NO.  OF  COUNTS 

CLR 

R 2 

RESAVE  R.2  TO  0 

OEC 

r: 

STORE  ONE  DATA  POINT 

RTI 

f 

DISMISS  INTERRUPT 

FINISH 

PRINT 

# M S G 4 

PRINT  OUT  CONTENTS  OF  MSG 4 

CLOSE 

#1 

CLOSE  CHANNEL  #1 

EXIT 

MSG1  . 

ASC  IZ 

/NO.  OF  BLOCH  C = Z5 

>:•  NO  OF  CYCLES'/ 

EVEN 

user 

. ASC  IZ 

/BEGIN  THE  DATA  AC OU I S I T I ON ' 1 / 

EVEN 

MSGS  . 

. ASC  I Z 

/STORE  €•  4 0 0 DATA  PC 

il  NTS.’ 

EVEN 

MSG  4 

. ASC  IZ 

/FINISH!'  / 

. EVEN 

BLOCK1 . 

. WORD 

0 

BL0CK2 . 

. WORD 

0 

CYCLE . 

WORD 

0 

LREC  : 

WORD 

0 

AREA 

BL  Ml 

5 

NAME 

RAD50 

/Ok  SHIH50DA1  ' 

DFiTH  STORED  IN  DISK  AS  THIS 

EVEN 

F I LE-NhME 

BUFFI 

. BL  KM 

14400 

BUFFI?: 

. BL F W 

14400 

HNOR  = 

CSECT 

TELETT 

DATA  IN 

BLEW 

40 

DA' A TRANSFER  FROM  ' T T I N' 

CSECT 

END 

MAIN 
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CSECT 

TITLE-  ASCIZ  COC-E  TO  DECIMAL  VALUE 


. GLOBL 

TTIN 

. MCALL 

. . V2.  . 

. REGDEF 

. . V 2 . . , . REGDEF 

TTIN 

MOV 

RO,  -<SF  > 

MOV 

R 1 . - ( S P > 

MOV 

R 2 . - < S P > 

MOV 

PS. -<SP> 

MOV 

R4, -CSP> 

MOV 

R5. -<SP> 

MOV 

#TEMF\ R2 

INLOOP  : 

. TTVIN 

•:  r 2 > * 

BICE: 

#200. RO 

CMP  6 

#15,  RO 

BNE 

I NLOOP 

MOV 

#DATA  I N,  R4 

CL  R 

F:1 

DELTA . 

CLR 

R 4 ) ♦ 

INC 

R 1 

C Mr1 

#20,  R 1 

BNE 

DELTA 

MOV 

WDATAIN- R4 

CLR 

RO 

MOV 

# T E M P . R2 

CMPB 

< R 2 > . #12 

BNE 

L OOP0 

INC 

F 2 

LOOPO 

CLR 

F 1 

LOO  PI 

CLR 

F 5 

CLR 

r: 

BICB 

# : 0 0 • ( F 2 

CMPB 

#15.  ( P 2 > 

BNE 

AGAIN 

MOV 

#100.  RO 

B F: 

LOOPS 

AGAIN 

CMPB 

#54, • R2> 

BED 

LOOPS 

CMP 

#0.  Ri 

BNE 

L 0 Ci  F 2 

BICB 

#S0.  < R 2 -• 

MOVE. 

( F;  2 > ♦ ■ FI 

BR 

L OOF'  1 

l OOF  2 

ADD 

Rl.  F5 

TTVIN 

; PUSH 

RO 

ON 

STACK' 

i PUSH 

Rl 

ON 

STACK 

. PUSH 

R2 

ON 

STACK 

. PUSH 

RS 

ON 

STACK 

; PUSH 

R4 

ON 

STACK 

. PUSH 

R5 

ON 

STACK 

. INITIAL  I 

ZE 

R2  WITH  ADDRESS 

. INPUT 

DATA 

FROM  DECWRITER 

, I N ASCIZ  CODE/ CLEAR  UPPER  PIT 
. 'CR'  ENDS  INPUTTING  DATA 

. INITIALIZE  R4  WITH  ADDRESS 
, P1=0 

i CLEAR  CONTENTS  OF  R4 
. CLEAR  20  WORDS 


, A LF'  LEFT? 


. TRANSLATE  ASCIZ  CODE  TO 
, DECIMAL  NO.  AND  STORE  RESULTS 
; IN  THE  LOCATION  OF  DATA  IN 


-CLEAR  UPPER  TWO  BITS 
, A ' CR'  FROM  DECWR I TER? 

, NO,  GET  NEXT  DATA 
, VES/  FINISH  INPUT 
/BRANCH  TO  LOOPS 
; A FROM  DEC  WRITER? 

• VES,  BRANCH  TO  LOOPS 
, NO/  R i * 0 '■ 

,N0. BRANCH  TO  LOOPS 
. NO. SUBSTPAC T DATA  BV  SO 
, MOVE  DATA  TO  R1 
•BRANCH  TO  L00P1 
. MULTI  PL  V BV  10 


C3  - NEW101 .MAC,  Conversion  of  ASCII  Codes  to  Numerical  Values 


APPENDIX  C 


in: 

RS 

CMP 

#12,  RS 

BNE 

LOOPS 

BICB 

#260,  <.  RS  > 

CLR 

Rl 

MQVB 

< R 2 > ♦ < Ri 

ADD 

RS.  Rl 

BR 

LOOF'l 

LOOPS : 

MOV 

Rl.  ( R4  > + 

CMP 

4100,  RO 

BED 

OUT 

INC 

RS 

JMP 

L OOPO 

OUT  : 

MOV 

( S P > ♦■ » R5 

MOV 

<SP >♦,  R 4 

MOV 

< S P > ♦ * RS 

MOV 

< SP  >♦. RS 

MOV 

<SP  • + , PI 

MOV 

c S F' + . R 0 

RTS 

PC 

TEMP 

. BLKW 
. EVEN 

40 

. CSECT 

TELETT 

DATA  IN : 

. BLkW 
. CSECT 

40 

. END 

TTIN 

INCREASE  RS  BV  1 
R 2 » 1 0 ? 

NO. KEEP  ADDING 
CLEAR  UPPER  FOUR  BITS 
CLEAR  R1 

MOVE  INPUT  DATA  TO  R1 
ADD  THE  PREVIOUS  DATA 
BRANCH  TO  LOOP! 

SAVE  RESULT 
FINISH? 

VES 

INPUT  NEXT  NUMERICAL  NO 

JMF  BACK  TO  LOOP0 

POP  STACK  TO  R5 

POP  STACK  TO  R4 

POP  STACK  TO  P3 

POP  STACK  TO  R2 

POP  STACK  TO  Ri 

POP  STACK  TO  RO 

RETURN 


, SPEC  I AL  LOCATION  FOR 
, RESERVED  INPUT  DATA 
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. C SECT 

TITLE  • WAVEFORM  SEGMENTATION  AND  RESTORftGE 


. GLOBL 

MAIN,  TTIN,  GRIDV, 

DISP,  CMPR1 , CMPR2, CMPR3 

GLOBL 

TIMEX, TIMEV 

MCALI. 

V2.  . , REGDEF, 

FETCH,  . ENTER,  . LOOKUP 

MCALL 

. V?.  . 

REGDEF 

TTVIN,  PRINT, . READW,  WRITW,  . CLOSE, . EXIT 

MULTI  A,  6 

, A = 2 5 6 *6 

. MACRO 

CL  R 

RO 

•EACH  BLOCK  CONTAINS  256 

CL  R 

A 

, WORDS 

ADD 

#400. A 

INC 

RO 

CMP 

B.  RO 

BNE 

. -14 

. ENDM 

MAIN 

. PRINT 

# M SGI 

•PRINT  OUT  CONTENTS  OF  MSG1 

JSR 

PC, TTIN 

.CALL  'TTIN'  FOR  DATA  INPUT 

MOV 

#DAT  A I N • P2 

, SET  UP  R2  WITH  ADDRESS 

MOV 

R 2 • *• , E’ i OC  F 1 

•MOVE  DATA  TO  BL0CK1 

FETCH 

# HNDR 1 ■ # NAME  1 

, DEFINE  FILE 

LOOKUP 

# A R F A 1 . #1 . # N A M E 1 

, LOOKUP  EXISTING  FILE 

. FETCH 

#HNDR2.  # N A M E 2 

, DEFINE  FILE 

ENTER 

# H R E A 2 ■ #?,  # N A M E 2 

, #-l  , ENTER  A NEW  FILE 

CL  R 

TIME! 

.SET  UP  TIMING  REFERENCE 

Cl  R 

TIME  2 

.EQUALS  ZERO 

SUB 

#4.  PL OCt  1 

• INITIALIZE  BL0C-K1 

AL  PHA 

ADI* 

#4.  E: L OCkl 

•TRANSFER  DATA  FROM  DISK 

REAf'W 

# H R E H 1 . #1  . # B l.i  F F S 

0,  #2000,  BL0CK1 

MOV 

#0.  F I PST 

, INITIALIZE  FOR  DISPLAV 

MOV 

#2000. LAST 

, INITIAL  IZE  FOR  DISPLAV 

BETA 

JSR 

PC, GRIDV 

-DISPLAV  C'AT A ON  SCOPE  BV 

JSR 

PC, DISP 

• FREOUENC  V VERSUS  TIME 

. PRINT 

# M S G 2 

•PRINT  MESSAGE 

JSR 

PC. TTIN 

.CALL  'TTIN'  FOR  DATA  INPUT 

MOV 

#[ AT AIN. R2 

, SE T UP  R2  WITH  ADDRESS 

MOV 

t R 2 > ♦ , PS 

•MOVE  INPUT  DATA  TO  RS 

CMP 

PS,  #1 

, R 3 = i 

PEG 

AL  PHA 

, VES. READ  NEXT  1024  DATA 

CMP 

#4.  RS 

, R S 24’’ 

ELE 

DELTA 

• VES, GOTO  DELTA 

MOV 

CSS"**.  FIRST 

i VES. INPUT  DATA  TO  FIRST 

MOV 

< F 2 *♦,  LAST 

-INPUT  DATA  TO  LAST 

CMP 

# 2 . R S 

. F:  S = 2 ? 

BEG 

BETA 

. VES, RE-DI SPLAV  SOME  DATA 

DELTA 

MOV 

#6LIF  F SO  ■ R4 

•SET  UP  R4  WITH  ADDRESS 

MOV 

F 4 , R 5 

, R5=R4 

ADD 

#4000, R5 

, F 5 = R 5 ♦ 2 0 4 6 

ASL 

F I C'ST 

, FIRST=FIFST42 

ADD 

FIRST,  R 4 

, R4= ADDRESS  OF  FIRST  DATA 
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GAMMA 

C MF 

• F 4 ’■  , 4-100 

E NE 

PHI 

ROO 

4 I , R 4 

MOV 

v R 4 •♦•  T I M E 4 

C M r' 

F 4 F 5 

E:GS 

EPSIL 

J 3 R 

FC,  T I M e : : 

PH  I 

MOV 

< R 4 • ♦ , T I M t 4 

J 3 R 

PC.  T ! M E V 

C HP 

F 4,  F 5 

PL  T 

GAMMA 

EPSIL 

C MF 

4 4.  F' : 

BNE 

THETA 

JMP 

C' ME  GA 

THETA 

PRINT 

4 MS  GO 

J S R 

P C > T T I N 

MOV 

4 i*  A CAIN.  PC' 

MOV 

• F , E l Or  1 :• 

CMP 

4 0.  BLOCK I 

5 E 0 

MO 

ROD 

4 4 . F L 0 . K 1 

MUL"I 

L F FC  0.  BLOC  1 

PEAOU 

4 AC  £ Hi  ■ 41  . lit  Ur 

SUB 

Hi.  [ i OC  1 1 

AOO 

EL  OCI  . BLOCK  ] 

MOV 

1 BUFFI. R4 

C l R 

F 0 

MORE  1 

CMP 

• F 4 • ■ 4-100 

PNE 

m ore: 

A 0 0 

4 . R 4 

H OV 

• F 4 4,  r I ME  4 

J SF' 

Ft , time:: 

Af'O 

4 j . F 0 

MORE  2 

MOV 

' R 4 . T I Mi  E 4 

J S P 

PC,  TIMES- 

I NC 

F 0 

C HP 

F 0,  l F Et  0 

EL  1 

HC  RE  1 

MU 

JMP 

AL  PHH 

, TIME 

REFERENCE  HAS  Ft 

OMFGA 

F R I N T 

4MSGZ 

PRINT 

4 MS  A 4 

J S R 

F'C  • T T I N 

MOV 

• OAT  A I N,  PC 

MOV 

'.  PI  '•  + ■ PL  OC  1 r 

MOV 

( PC  ■ ♦ . TOTAL 

MOV 

(FI  ♦•block: 

MOV 

• FI  •♦.  CS'CLEl 

MULTI 

LRECi.  B l 0 0 i < 

MUL  T I 

L R E L I E:  L 0 1 1 . 

MOV 

• BUFFI. R 4 

. t R4  ;•  * - 1 00 '' 

. N 0 , G 0 T C1  F'  H I 
.POINT  TO  NEXT  OATA 
, SAVE  OATA  IN  TIME4 
, FINISH'' 

■ V E S 

• PcAL-TINE  OF  UNITING  LOOP 
. SAVE  OATH  IN  TIME 4 
•REAL-TIME  OF  COUNTS 

• FINISH"' 

• NO.  GOTO  GAM M A 
, S' E S / RI  =4  "' 

NO 

. VES-REAO  NEXT  1014  C- ft T R 

, PRINT  MESSAGE 

-CELL  ' T T 1 N ' FOP  INPUT  OATA 

•SET  UP  PC  WITH  ADORE SS 

• MOVF  INPUT  TO  BLOCKS 

. PL  OC  f 0 = 0'' 

• VES 
, NO 

, L?ECI  = ISt"’  BLOCK  2 
-i  ■ IFECI. ELOCPi  / REAP  ORTA 
, FROM  01 SF 
•SET  UF  BL0CK1 
. SET  UP  R4  WITH  AOCRESS 
, R 0 = 0 

• lP4  ' = -100.  *■ 

• NO 

•VES. POINT  TO  NEXT 
/ S A v E C'  A T A IN  T I II E 4 
•REAL-TIME  OF  WAITING  LOOP 
/ F' H S S T NO  0 r i T A 
.SAVE  OATA  IN'  TIME4 
•REAL-TIME  OF  COUNTS 
.•  R 0 = R 0 ♦ 1 
/FINISH"' 

NO 

.VES- JUMP  PACK. 

EN  SET  AS  TIMEl  ANO  ' TIMES' 
FRINT  MESSAGE 
•PRINT  MESSAGE 

• CALI  1 T T I N ' FOR  OATA  INPUT 
•SET  UP  PC  WITH  AOORESS 
•MOVE  INPUT  TO  BLOCKI 
•MOVE  INFUT  TO  TOTAL 
•MOVE  INFUT  TO  BLOCKS 
•MOVE  INFUT  TO  CVCLE1 

• LREC 1 =00  S ♦ BLOCK  1 

• L RFCI*ISS ♦ BLOCKS 

• SET  UP  F 4 WITH  AOORESS 
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i TRANSFER  DATA  FROM  DISK  TO  MEMORV  WITHOUT  DISPLAV 
ROD  * #4.  BLOCKi  ; SET  UP  NEW  BLOCKI 

READW  4AREA1, #1- R4, LREC1. BLOCKI  .READ  DATA 

MOV  R4,  RESERi  .SET  UP  RESERi  WITH  R4 

MOV  R4,  RESER4  , SFR  RESER4  WITH  R4 

CLR  RESER2 

-THE  FOLLOWING  PORTION  OF  THIS  PROGRAM  IS  TO  STORE 
, ORTfl  6 S'  DEFINITE  NUMBER  OF  BLOCKS 
PRINT  4MSGS  .PRINT  MESSAGE 

CLR  BL0CK4  ;BL0CK4=0 

SUB  BLOCKS. BLOCK4  .INITIALIZE  BLOCK4 

TRV.  J S R PC.CMPR1  .FIND  BEGINNING  PART 

CMP  4200. FLAG  .END  OF  FILE? 

BED  FINISH  . VES 

J5R  PC.CMPP?  .FIND  CENTRAL  PART 

CMP  4200. FLAG  -END  OF  FILE? 

BED  FINISH  .VES 

JSR  PC.CMPP?  .FIND  BEGINNING  PART  OF  NEXT 

CMP  4200. FLAG  .END  OF  FILE? 

BED  FINISH  , VES 

ADD  Bl  OCRS,  BLOCK  4 .NO.  SET  UP  BLOC  K4 

i TRANSFER  DATA  TO  THE  DISK 
. WPITW  4 A R E A 2 . 42. 4BUFFS. L FEC2. BL0CK4 
PRINT  4 M S G 6 .PRINT  MESSAGE 

DEC  CVCLE1  i SAVE  ONE  WAVEFORM 

BEO  FINISH  .FINISH? 

J MP  TRV  ; NO. JUMP  BACK 

FINISH  PRINT  4MSG7  .FRINT  MESSAGE 

CLOSE  41  .CLOSE  CHANNEL  1 

CLOSE  42  .CLOSE  CHANNEL  2 

.EXIT  .FINISH1 

MS  AO  ASCIZ/HOW  M A N V BL  OCK  S(  = 224  > NC'T  TO  BE  DISPLAYED?/ 

. EVFN 

MSGi  ASCI  Z/WHAT  IS  1ST  BLOCK  OF  SHHSO.DAT  TO  BE  READ?/ 

EVEN 

MSG?  . . ASCIZ/SKIPd  • -RED  I SPL A V < 2 > -L AUNCH.  AFTER < S.  5>?0K<4>?/ 

. EVEN 

MSG 3 . . ASCIZ/START  DATA  PROCESSING'/ 

. EVEN 

MSG 4 : . ASCIZ/BL  OCf  ?■  =24,*  - TOTAL  ? BLOCKSOD?  CVCLE1  ?/ 

EVEN 

MSG5  . MSCJZ/SAVE  COMPLETE  WAVEFORMS  FOR  PROCESSING!/ 

EVEN 

MSG 6 . ASC IZ/TRANSFER  ONE  CVCLE  OF  DATA  TO  DISK!/ 

. EVEN 

MSG?  . ASCir/FINISH*  '/ 

. EVEN 

NAME  1 . FAD50/DK  SHI H50DAT.'  .WHERE  DATA  STORED  BV  100KHZ 

. Oc  REAL-TIME  CLOCK 

NAME?  PAD50/DK  SHIHlOrfiT/  .DATA  SAVED  FOR  PROCESSING 
CVCLE1  WOF  D 0 
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WORD 

'0 

BLOCf  4 

. WORD 

o 

L R E C 2 

WORE 

0 

AREF i2 

Eli  w 

c 

HNPR1 

E'LI  W 

5 00 

HNDR2 . 

.BLEU 

500 

. c sec  r 

TSLETT  SECTION 

FOR  A SC  12 

DAT  ft  in 

. BL  L IJ 

40  -CODE  TO 

DECIMAL 

C SEC  T 

V DATA  . SEC  T ION 

FOR  Y-AXIS 

V F 051 

BLEW 

VP  05  2 

E LL  W 

“ 

vr  ns: 

B L 1 W 

_ 

VS CALF 

WOP  D 

0 

CSEC  T 

HD  AT  A section 

F 0 R X - A X I S 

X F ij  3 1 

E:  L t W 

... 

X FO  S 2 

BL  1.  W 

XF  032 

E:  L 1 W 

X SCALE 

WORD 

TEST 

WORD 

CSECT 

NUMBER  -SECTION 

FOR  NUMERICAL 

NO 

E-ATE 

76. 111- 111- 105  76 

CHARACTERS 

re 

BYTE 

0-102  1 7 7 , 1 0 0 . 0 

n: 

. E VTF 

142-  1 ; l.  1 1 - 105-  102 

n: 

BYTE 

4 2,101  ill.  lll.ee- 

N4 

BYTE 

2 0,  24,  22,  1 77.  20 

N5 

BYTE 

4 T i Ci  5 • 1 0 5 , 1 0 5 . 7 1 

NO 

BYTE 

76.  111.  Ill,  111, 62 

NT 

BYTE 

101 . 41. 21. 11, 7 

NS  . 

E-  V T F 

66. Ill, 111, 111. 66 

NS* 

BYTE 

4c'.  1 11,  111,  111,  76 

SPACE  . 

BYTE 

0.  0.  0.  0,  0 

EVEN 

C SEC  T 

S C 0 F'c  -SECTION 

FOR  DI5RLAV 

E:  IJ  F P 3 0 

E LL  u 

2 0 0 0 

FIRST 

WORD 

LAST 

WORD 

C SEC  T 

SECOND  -SECTION 

FOR  TIMING 

TIME1 

WORD 

0 

TIME 2 . 

WORD 

0 

Tint  4 

WORD 

Ci 

CSEC  T 

TCNFR  SECTION 

F 0 R W h V E F 0 R M 

DAT A NO 

WORD 

0 SEGMENTATION  AND  STORA 

Pi  R E A 1 

ELL  W 

c 

FLAG 

WORD 

0 

l.PECl 

WORD 

0 

TC  TAL . 

WORD 

0 

BLOC-  1 

WORD 

o 

BLOCK  2 

HOP  D 

1.1 

F c 5 E P 1 

WORD 

c 

PE  SEP  2 

MORI 

0 

Pt ser: 

WOP  D 

0 
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RESER4.  WORD 

buffs  .blkw  10 

BUFFI:  BLKW  i4OO0 

BUFF  3 : .BLKW  3500 

. CSECT 

. END  MB  IN 
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. CSECT 
. TITLE 

N E W 250.  MAC 

. GLOBL 

TIMEX 

. M CALL 

. . VS.  . , . REGDEF 

TIMEX . 

. . V u 

. REGDEF 
MOV 

RD, -< 5P> 

MOV 

Ri, -<SP> 

MOV 

PS, -<SP> 

MOV 

F:3,  -(SP) 

MOV 

R4 , -CSP) 

MOV 

R5,  - ( S P > 

NEC 

TIME  4 

MOV 

TIME 4, R4 

CLR 

Ri 

CL  R 

RS 

CMP 

#0,  R4 

BLT 

ALPHA 

JMP 

OUT 

ALPHA : 

CMP 

#50  , P4 

BLT 

MORE  1 

BETA 

ADD 

#655  , RI 

DEC 

R4 

CMP 

#0,  R4 

BLT 

BETA 

MOV 

T I ME 4 . R4 

GAMMA 

SUB 

#1000.  , Ri 

INC 

TIMEi 

CMP 

#1000  , RI 

; (TIMED.  < TIMES >=Ev  65535* 
, H B S < T I M E 4 > 


, PUSH  RO 
; PUSH  Pi 
, PUSH  RS 

, push  r:< 

, PUSH  R4 

, PUSH  R5 

. T IME4  = -T I ME4 

, SET  UP  P4  AS  T I ME 4 

. R 1 = O 

; RS*0 

; 0 < R 4 ? 

> VES. 

; NO, JUMP  TO  RETURN 
, 50CR4? 
i VES. 

, R1=655*T I ME4 
. R 4 = R 4 - 1 
, R4<0? 

, VES, 

; SET  UP  R4 
, TIMF1*R 1/100O 
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BLE 

GfiMMfl 

JMP 

ETR 

MORE1 : 

RDD 

#€55.  , R1 

ROC 

R2 

DEC 

R4 

CMP 

#0,  R4 

BLT 

M0RE1 

M0RE2: 

SUB 

#1000.  , Rl 

SBC 

R2 

INC 

TIME1 

CMP 

#0.  R2 

BLT 

MORE  2 

CMP 

#1000.  . Rl 

BLE 

M0RE2 

ETR: 

MOV 

TIME4, R4 

RDD 

Rl/ TIME 2 

CLR 

Rl 

MORES : 

RDD 

#35.  , Rl 

DEC 

R4 

CMP 

#0.  R4 

BLT 

MORES 

M0RE4 : 

SUB 

#1000.  . Rl 

INC 

TIME2 

CMP 

#500.  / Rl 

BLE 

MORE  4 

OUT: 

MOV 

<SP>*. R5 

MOV 

<SP>+/  R4 

MOV 

<SP>*z  PS 

MOV 

<SP>«-z  R 2 

MOV 

(SP)*/ Rl 

MOV 

<SP>*,  RO 

RTS 

PC 

. CSECT 

SECOND 

TIME1 : 

. WORD 

TIME2: 

. WORD 

TIME4: 

. WORD 
. CSECT 

. END 

TIMEX 

; ( R2  > . <Ri>»655*FiBS<TIME4> 
/ DOUBLE-PRECISION 


; T I ME1 * ( R2  > . (RIVIBOO 
; DOUBLE-PRECISION 


i SET  UP  R4 
i SAVE  R1  TO  TIME2 

i Ri“35*RBS(TIME4> 


i TIME2*(R1>/1000 


, POP  R5 
, POP  R4 
, POP  R3 
; POP  R2 
i POP  R1 
i POP  R0 
; RETURN 

> SECTION  FOR  TIMING 
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. CSECT 

. TITLE 

TIME  REFERENCE  OF  lOOK  HZ 

. GLOBL 

TIMEV 

. MCALL 
. . V2.  . 

. REGDEF 

. . V2.  . . . REGDEF 

TIMEV : 

MOV 

R0> -<SP> 

PUSH  RO 

MOV 

Ri/ -<SP> 

PUSH  Rl 

MOV 

R2, -<5P> 

PUSH  R2 

MOV 

R?.  -<SF) 

PUSH  R3 

MOV 

R4,-<SP> 

PUSH  R4 

MOV 

R5, -<SP> 

PUSH  R5 

MOV 

TIME4, R5 

MOVE  DATA  TO  R5 

ALPHA : 

CL  R 

RO 

RO=0 

CMP 

#0,  R5 

R5<=6? 

BLE 

BETA 

VES, 

MOV 

#328.  , RO 

NO- SET  UP  RO 

BIC 

#100000- R5 

CLEAR  CARRV  BIT 

BETA: 

CMP 

#50.  , R5 

50>R5? 

BGT 

PHI 

VES. 

GAMMA: 

SUB 

#100.  . R5 

NO, RO=R5/10O 

INC 

RO 

CMP 

#50.  , R5 

BLE 

GAMMA 

DELTA. 

ADD 

RO- TIME2 

ADD  RO  TO  TIME2 

CMP 

#1000.  - TIME2 

10OOT1NE2? 

BGT 

PHI 

VES,  FINISH. 

NU: 

SUB 

#1000.  - TIME2 

NO, T I ME1=T I ME2/1000 

INC 

TIME1 

CMP 

#1000.  - TIME2 

BLE 

NU 

TIME=TIME1+TIME2/1000 

PHI  : 

MOV 

<SP>+-  R'5 

POP  R5 

MOV 

<SP>+, R4 

POP  R4 

MOV 

( S P ) + - R 3 

POP  R3 

MOV 

<SP>  + , R2 

POP  R2 

MOV 

<SP>+, Rl 

POP  Rl 

MOV 

<SP>+. RO 

POP  RO 

RTS 

PC 

RETURN 

. CSECT 

SECOND 

SECTION  FOR  TIMING 

TIME1  : 

. WORD 

TIME2 : 

. WORD 

TIME4 : 

. WORD 
. CSECT 

. END 

TIMEV 

C6  - NEW240.MAC,  Determination  of  Real  Time  from  Data 


APPENDIX  C 


73 


. C S E C T 

. TITLE 

BEGINNING  PORTION 

. GLOBL 

C-MPRl . TIMES',  TIME 

. MCfiLL 

. . VC.  . , . REGDEF,  . 

READW,  . PRINT 

. . v cT . . 

. REGDEF 

MACRO 

NX BUFF 

, MACRO  CALL  FOR  READING 

. PRINT 

8'ISGl 

DATA  FROM  DISK 

MOV 

8BLFF1, R4 

ADD 

BL0CK2,  BL  0CK1 

CMP 

TOTAL,  BL0CK1 

; TOTAOBLOCKi? 

BGE 

. + t? 

, VES, END  OF  FILE 

JMP 

FIN1 

. RE  RDM 

8AFEA1 , #1,  R4,  LREC1,  BL  0CK1  , READ 

CL  R 

p: 

, RE  = 0 

. ENDM 

CMPR1  : 

MOV 

F:0,  -<  SP) 

, PUSH  RO 

110  V 

FI.  -<sp) 

, PUSH  Ri 

MOV 

R2, -<SP> 

, PUSH  R2 

MOV 

RE, -CSP) 

, PUSH  RE 

MOV 

R4, -CSP) 

, PUSH  R4 

MOV 

F’5,  -CSP) 

; PUSH  R5 

MOV 

RESER1, R4 

; SET  UP  R4 

MOV 

PESER2, RE 

; SET  UP  RE 

MOV 

812, CATANO 

> SAVE  iO  F’TS  FOR  TIMING 

CLR 

TESTE 

; T E S T E = O 

ALPHA  : 

CL  R 

TEST1 

; TE  ST  1 = 0 

CLR 

TEST2 

; T E S T 2 = 0 

BET  A: 

INC 

R E 

, NEXT  DATA  POINT 

CMP 

RE. LRECi 

, END  OF  BUFFER? 

BGT 

GAMMA 

; VES- GOTO  GAMMA 

JMP 

EPSIL 

i NO- GOTO  EPSIL 

g ft  n ri  a : 

MOV 

8BUFF5, R2 

, SE T UP  R2  WITH  ADDRESS 

SUB 

820, R4 

-POINT  TO  PREVIOUS  ADDRESS 

CLR 

R5 

-SAVE  8 DATA  POINTS  OF 

DELTA 

MOV 

C R 4 ) + , i.  k 2 ) *■ 

i PREVIOUS  BUFFER 

INC 

R5 

CMP 

8L0,  Rf> 

BN  E 

PELT  A 

NX BUFF 

-READ  DATA  FROM  DISK 

EFSIL: 

CMP 

81,  TESTE 

j TESTE*1? 

BED 

EPSILi 

. VES, GOTO  EPSILI 

CMP 

R4,  RESER4 

, R 4 = R E S E R 4 ? 

BNE 

zeta 

, NO, GOTO  ZETA 

MOV 

81,  TESTE 

, TESTS=1 

EPSIL1 . 

CMP 

C R 4 ) , 8-100 

; ( R 4 ) = 1 0 0 ? 

BNE 

OMEGA 

-NO,  GOTO  OMEGA 

MOV 

81,  TEST2 

, VES,  T E S T 2 c 1 

ADD 

82,  RE 

.PASS  TWO  DATA  PTS 

ADO 

82,  R4 

; POINT  TO  NEXT  PT 
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MOV 

< R 4 ) ♦ , TIME4 

JSR 

PC, TIMEX 

CMP 

RS. LREC1 

6LT 

OMEGA 

JMP 

GAMMA 

OMEGA 

MOV 

< R 4 ) , TIME4 

JSR 

PC,  TIMEV 

ZETA 

CMP 

40, <R4> 

BLE 

THETA 

ADD 

42,  R4 

CMP 

41, TEST2 

BEQ 

RHO 

MOV 

41, TEST2 

JMP 

BETA 

THETA : 

CMP 

41,  TEST2 

BEQ 

ETA 

CMP 

41250.  , c R 4 ) + 

BLE 

PS  I 

JMP 

ALPHA 

F'SI  : 

MOV 

41, TEST2 

JMP 

BETA 

ETA  . 

CMP 

4900.  , (R4>  + 

BLT 

RHO 

JMP 

ALPHA 

RHO  : 

INC 

TEST1 

CMP 

42, TEST1 

BEQ 

PHI 

JMP 

BETA 

PHI  : 

MOV 

4BUFFS, Ri 

CL  R 

<R1>  + 

MOV 

T I ME 1 , (Rl)t 

MOV 

T I ME 2,  (Rl  >♦ 

CL  R 

( R 1 ) + 

MOV 

R4 , RESER1 

MOV 

R3, RESER2 

CMP 

RZ:,  410 

BLT 

SIGMA 

SUB 

420, R4 

BR 

LAMBDA 

SIGMA 

CLP 

R5 

MOV 

4BUFF5, R2 

KAPPA  . 

MOV 

< R 2 > «■ , < R 1 ) «■ 

INC 

DATANO 

INC 

R5 

CMP 

410,  R5 

6 ME 

KAPPA 

ASL 

RS 

SUB 

PS,  R4 

LAMBDA 

MOV 

L R4  ) ♦»  ( R 1 1 ♦ 

INC 

DATANO 

CMP 

R4 , PESER1 

.•  SAVE  DATA  IN  TIME4 
-CALL  TIMEX 
, BUFFER  IS  FULL? 

; NO, GOTO  OMEGA 
, VES,  GOTO  GAMMA 
, GET  THE  REAL  TIME 
, CALL  TIMEV 
; 0<=(R4>? 

, VES, GOTO  THETA 
, POINT  TO  NEXT  FT 
; 1=TEST2? 

, VES, GOTO  RHO 
; TE5T2=i 
; GOTO  BETA 
, TEST 2 = 1? 

• V'  E S , G 0 T 0 ETA 
; FREQ.  CBO  HZ? 

; VES,  GOTO  F'SI 
; NO, GOTO  ALPHA 
, TEST2=1 
, GOTO  BETA 
, FREQ.  <111  HZ? 

, VES, GOTO  RHO 
i NO, GOTO  ALPHA 
, TEST1 =TEST1 +1 
, TFST1=2^ 

, VES, GOTO  PHI 
, NO, GOTO  BETA 
; SET  UP  Rl  WITH  ADDRESS 
, SAVE  ONE  0 
.STORE  TIMING  DATA 
.STORE  TIMING  DATA 
i SAVE  ONE  O 
, 5 A V E R 4 
, SAVE  r:< 

, RSCB? 

, VES,  GOTO  SIGMA 

> NO, R4  = R4 -16 
.GOTO  LAMBDA 
, R 5 = O 

> SET  UP  R2  WITH  ADDRESS 
i STORE  S DATA  PT5 

; IN  TEMPORAR V BUFFER 


• R S = F'  2 2 

. R 4 = R 4 - R S 

.STORE  DATA  POINTS 
, UP  TO  CURRENT  DATA  PT 
■WHICH  FITS  THE  CRITERION 
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BNE 

LAMBDA 

J I'i  P 

'-UT 

FIN1  : 

MOV 

#200, FLAG 

OUT 

MOV 

Rl, RESERS 

MOV 

( SP +,  PS 

MOV 

<SF 5 * , R 4 

MOV 

<sp  r: 

MOV 

(SP)-i-,  RE 

MOV 

< SP.-  + . Rl 

MOV 

(SP>+, RO 

RTS 

PC 

MSG1  . 

. ASCI2, 
. EVEN 

READ  AT  £/ 

TEST1 

. WORD 

0 

TEST2 

. WORD 

TESTS : 

. WORD 

0 

. CSECT 

TCMPR 

DATANO : 

. WORD 

0 

APEA1 

BLkll 

5 

FLAG 

. WORD 

0 

L REC 1 ; 

. WORD 

0 

TOTAL  . 

WORD 

BLOCK1 

WORD 

0 

BLOC  (2 

. WORD 

0 

RESEP1  . 

WORD 

0 

RESER2 . 

. NOR  D 

0 

RESERS : 

WORD 

0 

RESER4 . 

. WORD 

BUFFS : 

. BLEW 

10 

BUFFI  : 

. BLEW 

14000 

BUFFS . 

. BLKW 

S'  500 

. CSEC T 

SECOND 

T I ME  1 . 

WORD 

TIMET 

. WORD 

TIME  4 

. WORD 
. CSECT 

END 

CMPP1 

i LOOP 
, G0T0  OUT 
» Sf  T UP  FLAG 

, POP  R5 
, POP  R4 
, POP  R3 
, POP  R2 
. POP  R1 
< P 0 P R 0 
, RETURN 


> SECTION  FOR  WAVEFORM 
.SEGMENTATION  AMP  STORAGE 


, SECTION  FOR  TIMING 
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CSECT 

TITLE 

CENTRAL  PORTION 

. GLOBL 

C M P R 2 , TIMEV,  TIMEX 

. MCALL 

VJl 

. . V£  . . REGDEF. 

RE  ADM,  PRINT 

V «.  . . 

REG DEP 

. MACRO 

NX BUFF 

, MACRO  CALL  FOR  READING 

. PRINT 

• MS61 

, FROM  DISK 

MOV 

• BUFFI. R 4 

ADD 

BLOCK?, BLOCK! 

CMP 

TOTAL, BL0CK1 

, TOTAL>BLOCKi? 

BGE 

. *6 

. VES, FINISH 

J M P 

FIN1 

. REA  D W 

•AREA1 . *1 . R4, LREC1 ■ BLOCK!  , READ 

CL  R 

r; 

ENDM 

CMPR2 

MOV 

PO, - < SF  ) 

, POSH  RO 

MOV 

Rl. -v  SP  • 

i PUSH  Rl 

MOV 

R£> -<SP> 

• PUSH  R? 

MOV 

RS,  - < S P ' 

. PUSH  r: 

MOV 

R4 . - < S P > 

■ PUSH  R4 

MOV 

R 5 , - ( £ P • 

PUSH  RS 

MOV 

RE  SERI, R4 

S A " E R 4 

MOV 

RESERi'.  R S 

S A V E F S 

MOV 

R'  E ? E R S • R 1 

SAVE  FI 

MOV 

R : PC 

R Z - R I 

CLP 

P 0 

. R 0 = 0 

ALPHA 

CLR 

PS 

■ * f T 0 

6'cTA 

INC 

r: 

, f =r:  +i 

CMP 

r:. lfsci 

BUFFER  IS  FULL"' 

6GT 

GAMMA 

•VES. GOTO  GAMMA 

J M P 

E P 5 1 L 

, NO,  GOTO  EF'SIL 

GAMMA . 

MOV 

PFSEP2. R! 

. SAVE  R* 

MOV 

RESER1 • R4 

. SAVE  R4 

DELTA 

MOV 

< R 4 > + , t R i ) * 

•STORE  DATA  POINTS 

INC 

[ ' n T A N o 

INC 

r : 

. f ;■  = R S ♦ 1 

CMP 

r:, lrec i 

BNE 

DELTA 

. L OOP 

NXBUFF 

, READ  DATA  FROM  DISK 

EPS1L 

CMP 

< F 4 • #-100 

, (.  R 4 > » - 1 0 0 

ONE 

ETA 

, N 0 , G 0 T 0 ETA 

ADD 

*£,  RS 

, PASS  TWO  DATA  PTS 

ADD 

Dt,  R 4 

, GOTO  NEXT  DATA  F'T 

MOV 

< R 4 + . T ! ME  4 

, SA  'E  DATA  IN  TIME  4 

JSR 

Pi.  TIMEX 

• C». CL  TIMEX  FOR  TIMING 

ETA 

MOV 

> P4 > . T l ME  4 

, SAVE  DATA  IN  TIME4 

J SR 

PC  TIMEV 

.CALL  TIMEV  FOR  TIMING 

CMP 

• 1,  PO 

» R 0 = 1 '' 

BEO 

1 AF  F A 

VFS, GOTO  K AF  PA 

C IIP 

• i o : c , <.  r 4 •>  ♦ 

f v e * *:s>r  h:- 

BLE 

LAME  DA 

•VES, GOTO  LAMBDA 

JMP 

AL  PHn 

•NO. GOTO  ALPHA 

KAPPA 

CMP 

• i 0 T 0 . • F 4 ♦ 

, FREO  > A i"  HC ^ 
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BGE 

LAMBDA 

, V E S , GOTO  LAMBDA 

JMP 

ALPHA 

.NO, GOTO  ALPHA 

LAMBDA 

INC 

R5 

, R 5 = R 5 ♦ 1 

CMP 

#6,  R5 

, R 5 = f.  ? 

BEQ 

MU 

, VES. GOTO  MU 

JM  P 

BETA 

, NO, GOTO  BETA 

MU  . 

CMP 

#1,  RO 

. R 0 = 1 ? 

BEO 

NU 

, VES. GOTO  NU 

MOV 

#1.  RO 

, NO.  R O 5 1 

JMP 

ALPHA 

; GOTO  ALPHA 

HU  : 

MOV 

PS. RESER? 

, SAVE  RESER2 

MOV 

R4. RE SERI 

; SAVE  RESERl 

CMP 

RESER3, R1 

. RESER3*Ri? 

BNE 

TAU 

, NO, GOTO  TAU 

SUB 

R2,  R3 

. RS=R3-R2 

TAU  : 

ASL 

R 3 

; R 3 = R 3 * 2 

SUB 

R3,  R4 

SIGMA . 

MOV 

(.  R 4 , ( R 1 ) ♦ 

.STORE  DATA  POINTS 

INC 

DATANO 

CMP 

R4» RESERl 

, R 4* RESERl? 

BNE 

SIGMA 

, NO, LOOP 

JMP 

OUT 

i VES, GOTO  OUT 

FI  HI  : 

MOV 

#200. FLAG 

, SET  UF  FLAG 

OUT  . 

MOV 

PI,  RESER3 

; SAVE  R1 

MOV 

CSP7*-, R5 

, POP  R5 

MOV 

< S P ) ♦ , R4 

, POP  R4 

MOV 

CSP R3 

, P 0 P R 3 

MOV 

(SP)+, R2 

• POP  R2 

MOV 

CSP)  ♦,  R1 

, POP  Pi 

MOV 

CSPjt, RO 

POP  RO 

RTS 

PC 

, RETURN 

MSG1  : 

. ASCI" 

/READ  AT  C/ 

. EVEN 

. CSECT 

TCMPR 

-SECTION  FOP  WAVEFORM 

DAT AND 

WORD 

O 

, SEGMENTATION  AND  STORAGE 

AREA1 : 

. BLI  W 

5 

FLAG 

WORD 

0 

LREC1 : 

. WORD 

0 

TOTAL  . 

WORD 

BLOCKi  : 

WORD 

0 

BLOC  K2 

WORD 

0 

PESER1  : 

WORD 

0 

RESERC 

WORD 

0 

resers 

WORD 

0 

RESER4 

WORD 

BUFFS . 

. ELK W 

10 

BUFF 1 : 

BLkW 

14000 

BUFF3 . 

. BU  N 

3500 

CSECT 

SECOND 

, SECTION  FOR  TIMING 

TIME1  ; 

WORD 

TIME2 

. WORD 

TIME4  . 

WORD 

CSECT 

END  CMFR2 
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. CSECT 

.TITLE  FINAL  PORTION 

. GLOBL  CMPR3, TIMEV, TIMEX 

. MCALL  . . V2.  . , . REGDEF. . READU, . PRINT 

. . V2.  . 

. REGDEF 


. MACRO 

NXBUFF 

MACRO  CALL  FOR  READING 

. PRINT 

#MSG1 

DATA  FROM  DISK 

MOV 

#PUFF1. R4 

SET  UP  R4  WITH  ADDRESS 

ADO 

BLOCKS, BLQCK1 

CMP 

TOTAL. BLOCK1 

TOT  AL>BL0CK1? 

BGE 

. +c* 

VES,  END  OF  FILE 

JMP 

FIN1 

READ  W 

K A R E A 1 . #1, R4 . LRECi, BLOCK! 

CL  R 

R3 

. ENDM 

CMPR3 . 

MOV 

RO, -<SP> 

PUSH  RO 

MOV 

Rl, -<SP> 

PUSH  Rl 

MOV 

R2. -<SP> 

PUSH  R2 

MOV 

R3,  - ( s p :• 

PUSH  R3 

MOV 

R4, -CSP> 

PUSH  R4 

MOV 

R5, -LSP) 

PUSH  R5 

MOV 

RESERi, R4 

SAVE  R4 

MOV 

RESER2, R3 

SAVE  R3 

MOV 

RESER3, Rl 

SAVE  Rl 

MOV 

R3-  R.5 

R5  = R3 

ALPHA : 

CLR 

TE5T1 

T E S T 1 = 0 

CL  R 

TESTS 

TEST  2 = 0 

BETA: 

INC 

R3 

R 3 = R 3 + 1 

CMP 

R3,  L RE Cl 

BUFFER  IS  FULL? 

BGT 

GAMMA 

VES, GOTO  GAMMA 

JMP 

EPSIL 

NO, GOTO  EPSIL 

GAMMA : 

MOV 

RESERS, PS 

SAVE  R2 

MOV 

RESERI, R4 

SAVE  R4 

DELTA : 

MOV 

< R 4 > , (Rl)  + 

STORE  DATA  POINTS 

INC 

DATANO 

INC 

R3 

R 3 = R 3 1 

CMP 

R3, LRECI 

R3* LRECi? 

BNE 

DELTA 

NO, GOTO  DELTA 

CLR 

R5 

NXBUFF 

VES. READ  DATA  FROM  DISK 

EPSIL 

CMP 

< R 4 > , #-100 

(R4)=-100'- 

BNE 

SETA 

NO, GOTO  SETA 

MOV 

#1, TESTS 

T E S T 2 = 1 

ADD 

US,  R3 

PASS  TWO  DATA  PTS 

ADD 

#2,  R4 

POINT  TO  NEXT  DATA 

MOV 

< R 4 > ♦ , TIME4 

SAVE  DATA  IN  TIME4 

JSR 

PC. TIMEX 

CALL  TIMEX  FOR  TIMING 

ZETA 

MOV 

* R4), TIME  t 

SAVE  DATA  IN  TIME4 

JSR 

PC,  T IMEV 

CALL  TIMEV  FOR  TIMING 
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CMP 

#0,  < P4  > 

BLE 

THETA 

ADD 

# 2 # R 4 

CMP 

#1#  TEST? 

BE  Cl 

MU 

MOV 

#1#  TEST 2 

JMP 

B E T A 

THETA  : 

CMP 

#1. TEST2 

BEQ 

ETA 

CMP 

#1250.  - «.  R 4 > + 

BL  E 

XI 

JMP 

ALPHA 

XI  : 

MOV 

#1#  TEST2 

JMP 

BETA 

ETA  : 

CMP 

#900. , <R4>  + 

BLE 

MU 

JMP 

ALPHA 

MU  : 

INC 

TEST1 

CMP 

#2, TEST1 

BEQ 

NU 

JMP 

BETA 

HU  . 

MOV 

R2# RESER2 

MOV 

F: 4 » RESFR'l 

MOV 

R4# RESER4 

SUB 

#2. RESER4 

CMP 

RE  SERB#  F:  J 

BNE 

T AU 

S U B 

R5,  RI- 

TAU  : 

ASL 

RS 

SUB 

R 2 # R 4 

OMEGA ; 

MOV 

< R 4 > + # < R 1 > ♦ 

INC 

R5 

INC 

DATA  N 0 

CMP 

R 4 # PESERl 

BNE 

OMEGA 

BR 

OUT 

FIN1  : 

MOV 

#200# FLAG 

OUT  : 

CL  R 

C P 1 > + 

MOV 

TIME1# <R1>+ 

MOV 

TIME2#  < Rl>  + 

MOV 

#BUFF2# R1 

MOV 

DATflNO#  <.  R 1 ♦ 

MOV 

RESER2  R2 

CMP 

#10,  R2 

BLE 

PHI 

SUB 

R2» RESER2 

ASL 

R2 

SUB 

R2,  RESERl 

BR 

PS  I 

PHI  ; 

MOV 

R5» PE  SEP 2 

MOV 

R4 . RESERl 

, 0 = <(.,?  4 >? 

# VES# GOTO  THETFl 
/PASS  ONE  DATA  POINT 
; TEST 2=1? 

; VES#  GOTO  MU 
; TEST2-1 
; GOTO  BETA 
# TEST 2*1? 

# VES#  GOTO  ETA 

■ FREQ.  <90  HZ? 

# VES# GOTO  XI 

; NO#  GOTO  ALPHA 
# TEST2*1 
# GOTO  BETA 
; FREQ.  <110  HZ? 

# VES. GOTO  MU 

# N 0 # G 0 T 0 ALPHA 
# TEST1=TE5T1+1 

# T E S T 1 = 2 ? 

» V E S # G 0 T 0 N U 
# NO# GOTO  BETA 

# SAVE  RESER2 
; SAVE  RESERl 

# SAVE  RESER4 

# RESER4=RESER4-2 

# refers  =ri? 

# NO# GOTO  TAU 

■ R 2 = R 2 - R 5 

# R S * F.'  S ♦ 2 

i R 4 = F:  4 - R 2 

# STORE  DATA  POINTS 


; R4=RESER1? 

, NO#  LOOP 
; VES. GOTO  OUT 

# SET  UP  FLAG 
i SAVE  ONE  O 

, SAVE  TIMING  DATA 
# SAVE  TIMING  DATA 
-SET  UP  R1  WITH  ADDRESS 

# STORE  NO.  OF  DATA  PTS 

# R I = R E S E R 2 

# 3 < = R 2 ? 

# V E S # G 0 T 0 P H I 

# RE3ER2=RESER2-R2 

# R ^ = R 2 2 

, RESERl =RESER1-R2 
i GOTO  PS  I 
, SAVE  RESER2 
, SAVE  RESERl 
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SUB 

#10,  RESER2 

SUB 

#?0, RE SERI 

PSI  : 

MOV 

(SP>*,  R5 

MOV 

<SP>+, R4 

MOV 

< S P ) + , R 3 

MOV 

<SP>+, R2 

MOV 

<SF>+, R1 

MOV 

(SP)+, R0 

RTS 

PC 

HSG1  : 

ASCIZ/READ  AT  E/ 

. EVEN 

TEST 1 : 

. WORD 

TEST?: 

. WORD  O 

. CSECT 

TCMPR 

DATfiNO 

. WORD 

0 

AREA1  : 

. BLKH 

5 

FLAG: 

. WORD 

0 

LREC1 : 

. WORD 

D 

TOTAL  : 

. WORD 

BL0CK1 

. WORD 

0 

BLOCK? 

. WORD 

O 

RESER1 

. WORD 

0 

RESER2 

. WORD 

0 

RESER3 

. WORD 

0 

RESER4 

. WORD 

BUFFS : 

. BLKW 

10 

BUFFI  : 

. BLKW 

14000 

BUFF3  : 

. BLKW 

5500 

. CSECT 

SECOND 

TIHE1 : 

. WORD 

TIME?: 

. WORD 

TIME4  : 

. WORD 

. CSECT 

. END 

CM  PR  5 

; SET 

UP  RESER2 

, SET 

UP  RESER'l 

, POP 

R5 

; POP 

R4 

, POP 

RI< 

; POP 

R2 

; POP 

R1 

, POP 

R0 

, RETURN 

i SECTION  FOR  WAVEFORM 
; SEGMENTATION  AND  STORAGE 


. SECTION  FOR  TIMING 


09  - NEW440.MAC  (Continued) 


. C SECT 

.TITLE  GR I DV,  VSHOW  ;DISPLAV  GRIDS  OF  V-flXIS 

GLOEL  GRIDV. VSHOW  ; AS  FREQ. FROM  0 TO  200  HZ 

L F'  S V C = 1 T O 4 1 1* 

. MCALL  . . V2.  . / . REGDEF 
. V2.  . 

. REGDEF 
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. MACRO  SHOW  A, B,  C, 0 

, MACRO  CALL  FOR 

DISPLAYING 

MOV 

#A,  YF'OSl 

, ORTA  IN 

SUBROUTINE  'VSHOW 

MOV 

#0, VP0S2 

MOV 

»C,  VP0S3 

MOV 

#0,  VSCALE 

JSR 

PC , VSHOW 

. ENDM 

GRIPV . 

MOV 

RO, -<SP> 

, PUSH  RO 

MOV 

Rl, -<SP> 

, PUSH  Rl 

MOV 

R2> -<SP> 

, PUSH  R2 

MOV 

RS> -CSP> 

; PUSH  R3 

MOV 

R4. - < S F > 

, PUSH  R4 

MOV 

R5, -( SP> 

, PUSH  R5 

MOV 

#10000, LPSVC 

, ERASE  THE  SCOPE 

SHOW 

N2,  NO,  NO,  200 

i PISPLAV 

200  ON 

SCOPE 

SHOW 

Nl,  NS,  NO.  ISO. 

, PISPLAV 

ISO  ON 

SCOPE 

SHOW 

Nl,  NS,  NO,  ISO 

, PISPLAV 

160  ON 

SCOPE 

SHOW 

Nt,  N4 , NO,  140. 

, PISPLAV 

140  ON 

SCOPE 

SHOW 

Nl,  N2.  NO,  120 

, PISPLAV 

120  ON 

SCOPE 

SHOW 

HI,  HO,  NO,  1 00. 

, PISPLAV 

100  ON 

SCOPE 

SHOW 

SPACE,  NS,  NO.  SO. 

. PISPLAV 

SO  OH 

SCOPE 

SHOW 

SPACE,  Nt',  HO,  60. 

. PISPLAV 

60  ON 

SCOPE 

MOV 

<SP>*. R5 

, POP  R5 

MOV 

<SP>+, R 4 

, POP  R4 

MOV 

< S P > ♦ , P3 

; POP  R3 

MOV 

< S P > «■ , R2 

; POP  R2 

MOV 

<SP>+, Rl 

, POP  Rl 

MOV 

< SP  > +,  RO 

; POP  RO 

RTS 

PC 

, RETURN 

. CSECT  VPATA  > SECTION  FOR  V-RX I S 


VPOSi 

. 6LKW 

YP0S2 . 

. SL  kW 

j. 

YP03.' 

. BLKW 

VSCAl E 

. WORP 

0 

C SECT 

NUMBER  , SECTION 

FOR  NUMERICAL 

NO 

BYTE 

26,  121,  111,1  05-  ?6 

CHARACTERS 

Nl 

. SVTE 

0,  102,  127,  1 00,  0 

H2 

. BYTE 

142,  121,  1 11.  105,  102 

NS  : 

. SVTE 

42. 101, 111. 111. 66 

N4 

BYTE 

0 ■ 2 4 > 0 2,  1 . * » 2 0 

N5 

BYTE 

4 7,  105,  105,  105-  71 

N6 

BYTE 

76. Ill, 11 1, 111, 62 

NT  . 

BYTE 

101 . 41 , 21 , 11 • 7 

NS 

SVTE 

66, 111- ill- Ill- 66 

NO 

BYTE 

46. 1 1 1 • 111, 111, 76 

SFRCF 

BYTE 

EVEN 

. c sec r 

ENP 

0,  0 0 O-  0 

OF  I I V 
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. C SECT 

. TITLE  VSHOW 
. CL  DEL  VSHOW 
LPSVC=ir041S 
LPSvc;i=iro42e> 

LPSVCV=1 70422 
. MCALL  . V2.  . , . REGDEF 

. . V2.  . 

. REGDEF 
> 

.THIS  SUBROUTINE  DISPLAVS  NUMERICAL  CHARACTERS 
.ON  THE  SCOPE  AS  THE  INDEX  OF  V-AXIS 


VSHOW : 

MOV 

RO, -<SP> 

MOV 

Rl,  -<SP.> 

MOV 

R2, -CSP> 

MOV 

RS,  -<SP> 

MOV 

R4, -vSP> 

MOV 

RS. -<SP> 

CLR 

R 4 

CL  R 

TEST 

CLR 

L PS VOX 

CLR 

RS 

CL  R 

Rl 

ALPHA : 

ADD 

VSCALE, RS 

INC 

Rl 

CMP 

#2  0,  Rl 

BNE 

ALPHA 

MOV 

RS. VL INE 

SUE: 

#4  0,  P5 

MOV 

RS, VSTAR 

MOV 

VPOS1, RO 

f3R 

DELTA 

BETA 

MOV 

V P 0 S 2 , RO 

MOV 

V£  TAR,  PS 

BR 

DELTA 

GAMMA : 

MOV 

VF  OSS, RO 

MOV 

VSTAR,  R5 

DELTA . 

MOV 

#-5 . Rl 

EPSIL  . 

ADD 

#15,  R4 

MOV 

VSTAR, RS 

MOV 

#-7,  R2 

MOVE. 

CRO'1  ♦,  R 3 

ZETA  . 

ROLE: 

RS 

E:PL 

IOTA 

MOV 

#2002.  LPSVC 

ETA 

TSTB 

LPSVC 

PPL 

El  A 

MOV 

RS.  LPSVC  V 

MOV 

R4.  LF SVC X 

. PUSH  RO 
> PUSH  R1 
; FUSH  R'2 
. PUSH  RS 
, PUSH  R4 
; PUSH  RS 
. R 4 = 0 

; INITIALIZE  TEST 
.VALUE  OF  HONRI  =0 
; R 5 - O 
. R 1 = O 

. R 5 = F R E 6 * 1 6 


. NO.  LOOP 
. SAVE  FREC'*16 
. R 5 » R 5 - 3 2 
..  SAVE  RS 

.MOVE  FIRST  CHARACTER 
; GOTO  DELTA 

i MOVE  SECOND  CHARACTER 
; MOVE  HONRI.  GRID 
; GOTO  DELTA 
MOVE  THIRD  CHARACTER 
.MOVE  HONRI.  GRID 
. Ri*-5 
, R 4 = R 4 1 3 

.MOVE  CHARACTER  TO  RS 
. R 2 = - 7 

•MOVE  E PITS  TO  RS 
, RC'TAiE  0 El T TO  CARRY 
.CARRY  E;  I T = O ? 

> NO- SET  UP  STATUS  OF  SCOPE 
, SCOPE  READY? 

, NO.  WAIT 

. VES.  PUT  ONE  DOT  ON  SCOFE 
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IOTA . 

ft  PC’ 

#11,  R5 

, I NCREPSE  VERT  I.  BV  ? 

INC 

R2 

, TEST  7 DOTS  FOR  ft  COLUMN 

BNE 

SET  ft 

INC 

R1 

;5  COLUMNS  FOR  ft  CHARACTER 

BNE 

EPSIL 

INC 

TEST 

> TEST  THREE  CHARACTERS 

CMP 

#1,  TEST 

i Tt  ST  = 1 ? 

Bt  Q 

BET  ft 

iVES, TRY  SECOND  CHARACTER 

CMP 

#2, TEST 

i NO-  T E S T “ 2 ^ 

B 1 0 

GftMMft 

, YES, TRY  THIRD  CHARACTER 

ADD 

1*100,  P 4 

, HONRI.  INCREASE  BY  64 

MOV 

YLINE, R2 

KAPPA : 

MOV 

#2002, 1 PSVC 

-SET  UP  STATUS  OF  SCOPE 

MU  . 

TSTP 

L PS  VC 

, SCOPE  READY? 

BPL 

MU 

i NO, WAIT  FOR  READY 

INC 

F:4 

; YES, INCREASE  HONRI.  BV  1 

MOV 

R4,  L PSVC  X 

-GENERATE  ONE  DOT  ON  SCOPE 

MOV 

R2. LPSVCV 

-KEEP  SAME  VALUE  OF  VERT  I . 

CMP 

t*  r ? r t- , Fi  4 

, OUT  OF  RANGE  OF  SCOPE? 

BNE 

KAPPA 

, NO- KEEP  DISPLAYING 

MOV 

( SP> R 5 

• POP  R5 

MOV 

<SP>*,  R4 

, POP  R4 

MOV 

(SP  >♦,  RS 

, POP  RS 

MOV 

< S P > ♦ , R 2 

, POP  R2 

MOV 

<:sp->  + . ri 

, POP  Rl 

MOV 

C S P ) ♦ , R Ci 

, F'  0 F R 0 

RTS 

PC 

, RETURN 

YSTRR ; 

. WORD 

0 

YLINE  : 

. WORD 

0 

TEST 

WORD 

0 

. CSECT 

V D ft  T ft 

i SECTION  FOR  Y-AXIS 

VPOS1  . 

. BIKW 

\ 

VP0S2 

. Blk'M 

VPOSS 

. BLEW 

VSCflLF . 

. WORD 

0 

. CSECT 

. END  VSHOW 
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DISP  : 


DEVI 


START 
CRT1  . 


ALPHA 


CRT2 

READY 


. CSECT 
. TITLE 

DISP  THE  DATA 

OF  100k  H2 

. GLOBL 

DISP. XSHOW 

L PSVC  = 1 

70416 

LP5VCX* 

170420 

L P S V C V = 

170422 

. MCALL 
V2. 

. REGDEF 

. V2  ...  REGDEF 

MACRO 

SHOW  A.  B.  C 

i MACRO  CALL  FOR 

MOV 

#A.  XFOSl 

.DISPLAYING  X-AXIS 

MOV 

#B. XP0S2 

MOV 

#C.  XP0S3 

JSR 

PC.  XSHOW 

. ENDM 
MOV 

RO. -<SF> 

; PUSH  RO 

MOV 

R 1 . - ( S P > 

. PUSH  Rl 

MOV 

R2. -<SP> 

. PUSH  R2 

MOV 

r s . -(spy 

, PUSH  PS 

MOV 

R4. - < SP  > 

, PUSH  R 4 

MOV 

R5.  - < S P > 

, PUSH  R5 

MOV 

#BUF  F S O.  RO 

.SET  UP  RO  WITH  ADDRESS 

MOV 

FIRST. Ri 

• SET  UP  Rl 

MOV 

LAST,  R2 

, SET  UP  R2 

SUB 

Rl,  R2 

, R 2 = R 2 - R 1 = D A T A POINTS 

MOV 

R 2. DATAl 

, SAVE  NO.  OF  DATA  F'TS 

MOV 

* lDOOO, R4 

, R4=4096 

CL  R 

R5 

, R5  = 0 

SUB 

R2.  R4 

; R5  = 40S*<f/R2 

INC 

R5 

CMF 

R2 . R4 

BLE 

DEVI 

ASL 

Rl 

i P 1 = 2 *•  R 1 

ADD 

PI.  RO 

. R0= ADDRESS  OF  FIRST  DATA 

MOV 

FIRST, Rl 

. R 1 = N 0 OF  FIRST  DATA 

DEC 

Rl 

CL  R 

l PS vex 

, L P S V C X = O 

CMP 

> RO',  #-100. 

. VALUE  — 100? 

BNE 

ALPHA 

.NO.  GOTO  ALPHA 

ADD 

# 4 , R 0 

; SHIP  TWO  DATA  PTS 

ADD 

#2,  Rl 

MOV 

#S?00  . TEST 

i DISPLAY  IN  SPECIAL  FORM 

JMP 

MOVE  1 

MOV 

C R O '<>  ♦ • r 4 

.SAVE  DATA  IN  R4 

INC 

Rl 

. R 1 = R 1 ♦ 1 

CMP 

Rl ■ LAST 

, R12LAST? 

BLT 

CRT2 

, YES, GOTO  CRT2 

JMP 

OUT 

, GOTO  OUT 

MOV 

#2002. IPS VC 

.SET  UP  THE  SCOPE 

TSTB 

IF  SVC 

. SCOPE  READY?  . 
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GRMMRO  : 


GRMMfil  : 

DEV2  : 


MOVE  J 


OKO  . 


Ok  1 . 


Ok  01 


Ok  11 


Ok  21 
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BF'L 

REROV 

, NO, WRIT 

CMP 

#480  , R4 

, 4 SO ON 0.  OF  COUNTS? 

BLT 

GRMMRO 

; VES. GOTO  GRMMRO 

MOV 

#S5G0  .TEST 

> MSPLRV  IN  SPEC  I RL  FORM 

BR 

M0VE1 

, GOTO  MOVE 1 

CMP 

#0,  R4 

, 0<R4"> 

BLT 

GRMMR1 

, VES,  GOTO  GRMMR1 

MOV 

#200  - TEST 

, TEST  = 200 

JMP 

MOVE  1 

i GOTO  M0VE1 

MOV 

#85000- R2 

. DOUBLE-PRECISION 

MOV 

#:  o, p: 

CLR 

TEST 

SUB 

R4,  R2 

; TEST  = 1 . BOO, OO0/R4 

SBC 

RS 

INC 

TEST 

CMP 

#0-  PS 

BNE 

D E V 2 

CMP 

R4.  R2 

BLE 

D E V 2 

ROD 

R 5 .■  L PS  VCX 

i MSPLRV  R DOT  ON  SCOPE 

MOV 

TEST..  LPSVC  V 

MOV 

L F‘  s V CX,  X S C R L E 

CMP 

r RTRI , #150. 

, DRTR  PTSM50? 

EG  T 

Ok  01 

, VES, GOTO  0R01 

CLR 

RS 

MOV 

R 1 R 4 

SUB 

#12,  R 4 

, R 4 , ‘ 1 0 = P 2 

INC 

RS 

CMP 

#12. R 4 

BLE 

OCB 

CMP 

#0.  ?' 4 

, R 4 = 0 ? 

BEO 

Ok  1 

j VES,  GOTO  l1  k 1 

JMP 

CPT1 

, NO, GOTO  CRTI 

SUB 

#12,  F : 

, RS/10 

CMP 

#11.  R 3 

BLE 

Ok  1 

CMP 

# o , r : 

, R S * 0 ? 

BEO 

Ok  01 

, VES- GOTO  OKOl 

SHOW 

SPRC  E,  SPRC  E • S 

PRC  c .VERT  I.  LINE  ON  SCOPE 

JM? 

C FT  1 

GOTO  CRTI 

CMP 

#100  , PI 

BNE 

CM  1 1 

SHOW 

N 1 ■ NO. NO 

, MSPLRV  100  RND  VERT  I 

JMP 

CRT  1 

, LINE  ON  SCOPE 

CMP 

#200  PI 

BNF 

Ok  21 

SHOW 

N2 . NO. NO 

, DISFLRV  200  RND  VERT  I. 

JMP 

CRTi 

.LINE  ON  SCOPE 

CMP 

#‘00  , R 1 

BNE 

Ok : i 

SHOW 

N : . NO.  N 0 

.MSPLRV  :>00  RND  VERT  I 
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JMP 

CPT1 

i LIME  ON 

OK  SI  : 

CMP 

#400.  , R1 

BNE 

OK  41 

SHOW 

M4 » NO,  NO 

, DISFLAV 

JMP 

CRT1 

; L I ME  ON 

0K41  : 

CMP 

#500.  , PI 

SME 

OK.  51 

ShOW 

N5,  MO,  NO 

DISFLAV 

JMP 

CRT1 

, LINE  ON 

OF  51  : 

CMP 

# £ 0 0 , R1 

BME 

OK  61 

SHOW 

MB,  MO,  NO 

; DISFLAV 

JMP 

CRT1 

, LINE  ON 

OK  SI  . 

CMP 

#700  , R 1 

BME 

OK  71 

SHOW 

N 7 , MO,  NO 

, DISFLAV 

JMP 

CRT1 

■ LINE  ON 

OK  71 

CMP 

# 500  ,R1 

B N E 

of  ei 

S H 0 W 

MB,  MO,  MO 

; DISFLAV 

JMP 

C RT1 

, L INc  ON 

OK  Si 

CMP 

# ? 0 0 , R 1 

BME 

OF  Si 

SHOW 

MS,  NO,  NO 

; DISFLAV 

JMP 

C RT1 

LINE  ON 

OK  SI 

CM  P 

#1000  , R1 

B N E 

OK  10  1 

SHOW 

Ml , NO,  MO 

, DISFLAV 

OKI  Ot 

JMP 

C R T 1 

> LINE  ON 

OUT 

MOV 

( S P ..*  + > F 5 

; POP  R5 

MOV 

«'. SP>+. R4 

, F OP  R4 

MOV 

( i'  F > +• , R S’ 

, F u P'  R S 

MOV 

< SP  > ♦,  P i 

, POP  R£ 

MOV 

<SP>+,  P 1 

, FOP  R1 

MOV 

< S P • , R 0 

, P 0 P F.  0 

RTS 

PC 

, RETURN 

0 A T A 1 

WORD 

o 

CSECT 

SC  OF  s 

, SECT  I ON 

BUFFSO 

. BL'.W 

iC'OO 

FIRST  . 

. WORD 

LAST 

. WORO 

CSECT 

XDATrl 

, SECT! ON 

NF0S1 

BLKU 

S; 

KP0S2 

. El  KM 

7 

X P OSS  . 

. ELI  W 

XSCALE 

. WORO 

o 

TEST  . 

W 0 R 0 

CSECT 

NUMBER 

, SECTION 

MO  . 

. E'VTE 

7 B < 1 i 1 , 11 

1 , 105,  7 

B 

Ml 

E'VTE 

0 , 1 0 i 177 

. 100,  0 

Mi'  . 

6 V T £ 

1 4 i , 1 i 1 , 1 

1 1,  105, 

1 Oi 

C12  - NEL’330 .MAC  (Ccr.tlnued) 

SCOPE 


400  AND  VEPTI 
SCOPE 


500  AND  VEPTI 
SCOPE 


f.00  AMO  VEPTI 
SC  OPE 


700  AMO  VEPTI 
SCOPE 


COO  A NO  VEPTI 
SCOPE 


SOO  AND1  VEPTI 
SCOPE 


100  AMI1  VEPTI 
SCOPE 


FOR  DATA 


FOP  X-VALUE 


FOP  NUMERICAL 
l H A P A C T E R S 
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N3 

. BYTE 

N4 

. BYTE1 

N5 

. BYTE 

N6 

. BYTE 

N7 

BYTE 

NS 

BYTE 

N? 

. BYTE 

SPACE:  BYTE 

. EVEN 

CSfc'CT 

END 


42..  101,  111,111,  66 
3 0 , 24,  22,  1 7 7 , 20 
47,  105,  105,  105,  71 
76, 111, 111, 111, 62 
101,  41, 21,  11,  7 
66.  111,111,111,  66 
46,  111,  111,  111,  76 
0,  0,  0,  0,  0 


DISP 
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. CSECT 

. TITLE  K SHOW 
. GLOBL  XSHOW 
LPSVC=1704i6 
LPSVCX=1 70420 
L P S V C V = 1 7 0 4 2 2 
. MCALL  . . V2.  . , . REGDEF 
. . V2.  . 

. REGDEF 

; THI 5 SUBROUTINE  IS  TO  DISPLAY'  NUMBER  I CAL  CHARCTER* 
i ON  THE  SCOPE  AND  EXPAND  THE  HONRI.  SCALE 


XSHOH : 

MOV 

R0,  - ( SP  > 

; PUSH  RO 

MOV 

Rl, -<SP> 

, PUSH  Rl 

MOV 

R2, -<SP> 

, PUSH  R2 

MOV 

R 3 » - ( S P ) 

, PUSH  R3 

MOV 

R 4 , - ( SP) 

; PUSH  R4 

MOV 

R5, -<SP> 

; PUSH  R5 

CLP 

R4 

, R 4 = 0 

CLR 

L PS VCX 

, L P S V C X c O 

CL  R 

LPSVCV 

; L P S V C V = O 

MOV 

X 5 C ALE,  XL INE 

.SAVE  XSCALE 

SUB 

W26.  , XSCALE 

. XSCALE=XSC ALE-36 

MOV 

X p 0 S 1 , ^ 0 

.PIRST  CHARACTER 

BR 

GAMMA 

, GOTO  GAMMA 

ALPHA 

MOV 

X F 0 S 2 RO 

, 2ND  CHARACTER 

BR 

GAMMA 

, GOTO  GAMMA 

BETA  . 

MOV 

XROS 3, RO 

, 3RD  CHARACTER 

GAMMA 

MOV 

l»-5.  Rl 

, Rl =-5  FOR  5 COLUMNSXC HARAC 

DELTA: 

ADD 

♦♦14.  XSC  AL  E 

.POSITION  FOR  NEXT  COLUMN 

CLR 

R5 

, R 5 " 0 

C13  - NEW340.MAC,  Displaying  Horizontal  Coordinate 
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f 

I 


MOV  • 

# - 7 . R2 

MOVE: 

( F:  0 + , R 2 

SET  ft  . 

RC'LB 

R2 

BPL 

KAPPA 

MOV 

#2002, LPSVC 

PHI  : 

TSTB 

LPSVC 

BPL 

PHI 

MOV 

R5, L PSVCV 

MOV 

XSCALE, LPSVC 

KAPPA : 

ADD 

#il,  RS 

INC 

R2 

BNE 

ZETA 

INC 

R1 

BNE 

DELTA 

INC 

R4 

CMP 

# 1 , F:  4 

BEQ 

ALPHA 

CMP 

# 2 , K 4 

BED 

BETA 

MOV 

#120.  , F'5 

MU 

MOV 

#2002, LPSVC 

OMEGA: 

TSTB 

LPSVC 

BPL 

OMEGA 

INC 

R5 

MOV 

TES1 , R3 

SUB 

R5,  R3 

CMP 

#200, R2 

BLT 

PS  I 

ADD 

# £ 6 0 , R5 

MOV 

#40000, TEST 

PS  I 

MOV 

R5,  LF'SVCV 

MOV 

KLINE, L PS VCX 

CMP 

#7776,  F 5 

BGE 

MU 

MOV 

<SP>+, R5 

MOV 

<SP>  f,  F:4 

MOV 

<SP>+. F2 

MOV 

<SP>  + , p 2 

MOV 

< S F '>•*■,  F,  1 

MOV 

••  SP>-,  RO 

RTS 

PC 

XL INE . 

. WORD 

0 

. CSECT 

X DATA 

XPOSi  . 

. BLEW 

XF0S2 

BLKM 

XP0S3 . 

BLKM 

"SCALE . 

. WORD 

TEST 

. WORD 

CSECT 

. END  KSKC'iW 
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FOR'  7 ROWS/C  HfiRfiC. 
; SAVE  DPT  ft  IN  R2 
i R 0 7 ft  7 E R 2 

, CARRV  BIT  SET  , D I SPLAV  ft 
, SET  UP  SCOPE 
> SCOPE  RFADV? 

, NO,  WRIT 

, VES,  PUT  ft  DOT  OH  SCOPE 

, V-FOSITION  FOR  NEXT  DOT 
, R 2 * R 2 ♦ 1 

; FINISH  ft  COLUMN? 

; R1=R1+1 

, FINISH  ft  CHARACTER? 

, P 4 = R 4 + i 
, R 4 = 1 ? 

, V E S , TRV  2ND  CHARACTER 
, R 4 » 2 ? 

, VES,  TRV  3RD  CHARACTER. 

, R5-120 

..  SF  T UP  THE  SCOPE 
-SCOPE  READV? 

. NO.  WAIT 
; P5-R5+1 

, rl-»test*v-value 

, R2-R3-R5 
, 1 ft  2 R - ' 

, VES, GOTO  RSI 
, R 5 = R 5 + 3 8 4 
•SET  UP  TEST 
; VERTI.  STRAIGHT  LINE 

, 4 0 S 5 > = R 5 ? 
i VES,  GOT 0 MU 
, POP  F'5 
..  POP  P 4 
i POP  PS 
, F 0 F1  R 2 
, POP  R 1 
. POP  RO 
. RETURN 

i SECTION  FOR  X-VftLUE 
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THIS  PROGRAM  IS  TO  GET  THE  SLOPES  CORRESPONDING 
TO  THE  POSITIVE  AND  NEGATIVE  CONDUCT  I VI T I ES 
DIMENSION  MCI  792  >,  1 UNO < IDO),  TIMED1'  200  > , C0NDl20D> 
COMMON  M, G, TREF. IFRE01, IFREG2, IAXIS, ID1 
DATA  HAS  BEEN  STORED  IN  THE  DISK  UNDER  THE 
FILENAME  "SHIH10.  DAT" 

TVPE  :o 

FORMAT < IX, ' NO.  OF  BLOCK?. FIRST  AND  LAST  WAVEFORM  ?' ) 
ACCEPT  50,  I BLOCK,  I BEG  IN,  I FND 
FORMAT  * 5 I S > 

L R EC =256*  1 BLOCK 

IB  = IBLOCM  -.  IBEGIN-2) 

IFLAGZ:  = 0 
N«0 

DO  2000  J=IBEGIN, IEND 
IB=IB+IBLOCI 

SUBROUTINE  CHECK.  IS  TO  READ  DATA  FROM  DISK- 

CALL  CHECK  CM.  IB,  LREO 

ID1*5 

I D 2 = M < 1 >-€■ 

TVPE  100,  J,M(l) 

FORMAT  1 ////IK,  ' WAVEFORM  ' , 1 1 , 5X,  ' DAT A PTS  : ' , I 4 > 
SUBROUTINE  GRIDVtMAO  DISPI.AVS  V-AXIS  ON  SCOPE. 

CALL  GRIDS' 

SUBROUTINE  D ISP (MAO  DISPLAVS  DATA  POINTS  ON  SCOPE. 
CALL  DISPCM. I D 1 • ID2, IAXIS> 

TVPE  1 BO 

FORMAT  C1X,  EXPAND  ( 1 > '-'POS  CC-ND.  (2'-?NEG.  CGND.  (2>? 

N SK  I P ( 4 > ''STOP <.  5 > > 

ACCEPT  50,  I FLAG 

GOTO  a 90,  : 00,  400.  1 950,  2000  V.  I FLAG 
TVPE  200 

FORMAT < IX, ' LIMITS  OF  THE  X-AXIS  * V-AXIS  ?' > 

ACCEPT  50,  ID1,  I D 2 , IFRE01,  IFREQ2 
EXPANSION  OF  PART  OF  THE  WAVEFORM 
I F R E Q * = I M T C S 200.  /(IF  R E G 2 - I F P E C!  1 ; ) 

CALL  VEXFDC IFRE01,  IFREG2,  IFPEOD 
N 01  = 1 0 0 , O0O/FREG1, N02  = 100,  000/FRE02 
N 0 1 = I N T < < 1 0 O 0 0 . / I F R E Q 1 > + 10  +0  5> 

N 0 2 = I N T < a 0 0 0 O . / 1 F R E 0 2 S 1 0 tO  5 > 

I E X 1 = F R E C!  1 * 2 2 0 0 / < F R E 0 2 - F R E 0 1 > - 4 0 0 
IEX2-250, OOO/L  FREQ2-FRE01 ) 

I E X 1 = I N T ( I F R E 0 1 ♦•  < S 2 0 0 /(I  FRE02-  1 FRE01  > ) -400  > 

IEX2  = INT<  a 0000.  / cl  FRED2-I  FPEC1  > ->♦  10.  *2.  5> 

CALL  XEXFDOI,  IDi,  I D 2 , N 0 1 , N 0 2 • IEX1,  IEX2,  IAXIS> 

G 0 T 0 170 

IFvIFLAG:  EG  1 > GOTO  400 

tvpe  ;:.2-o 

FORMAT  ax,  ' M POSITIVE  C 0 N D I ! C 1 I V I T V # • ' ) 

SUBROUTINE  SLOF'l  IS  TO  FIND  POSITIVE  CONDUCT  I V I T V 
CALL  SL0P1 
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GOTO  600 
•100  TVF'E  420 

420  FORMATtlX.  ' ^NEGATIVE  CONDUCT  I VI  TV*#'  ) 

CALL  51  OF' 2 
IFLAGZ  *1 
COO  TVF'E  TOO 

?00  FORMAT <'  1 ' M SF  L flV  i\  CALC  <P0TH*1, ONLY  CALC 

» 

ACCEFT  50.  JFLAG 
GOTO  < IDO. 000. 1000>, JFLAG 
C R-.CRL>*5*E11.  RAM  I ARE  0225  FIND  5" 

C RESPECTIVELY  FOR  FLU NT  PROBE 

1000  G1 =1  24r*G«lE-i: 

TYPE  1100.  Gl.TRE" 

1100  F OR  MAT  i 4:-:.  'CONDUCTIVITY.  '.Eli*  4-16X. 'TIME  : ' . F 9 . 3 > 

N = N *-1 
I UNO  > N > «=  J 
TIM'EOv  N>*TREF 
con:-  N)«G1 
GOTO  170 

1550  TYRE  I960 

I960  FORMAT  v IX,  HON  MFlNY  WAVEFORMS  NOT  TO  EE  PROCESSED?') 
ACCE:'T  50.  IK 
I B - I B + 1 K *1  B L 0 C'  K. 

JaJ>  :k 
IFLAGI^O 
2000  CONTINUE 

3000  TYPE  3500 

3500  FORMA  '(1  OX.  ' WAVEFORM  . 5 X.  TIME',  SOX.  ' CONDUCTIVITY' 

<1  12, X.  (NO.'  . 6 X , • (SEO'  , 23X,  ' (MHO/CM)'  /’> 

TYPE  I 600.  ( I UNO  > I ). TIME 0 < I >, COND< I >•  I ° 1 • N > 

3600  FORMAT' 12K.  12,  5X,  F?  3,  1SX,  E12.  4.  ,‘> 

4000  STOP 

END 
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C HECK  . 


ftREA 
NANI . 

KNDR  = . 


. CSECT 

. TITLE 

'CHEC1 

GLOBL 

CHECK 

. MCftLL 

. . V2  . - PEGDEF 

. MCftLL 

. CLOSE. . EXIT 

VC.  . 

. REGDEF 

MOV 

PO  -cSP> 

MOV 

R 1 . - ( S F > 

MOV 

PC',  - ( 5 F'  > 

MOV 

PS. -<SP> 

MOV 

F.4,  - ( S P > 

MOV 

ft 5 - ( SP  > 

MOV 

2<P5>. PI 

MOV 

CM’.  R5>,  ft? 

MOV 

C*6(R5>.  r: 

FETCH 

ItHNDft.  * N ft M E 

LOOKUP 

It  ft  ft  E A • * 1 , It  N ft  M E 

REftDW 

it  ft  p e ft , #x,  pi,  p;. 

. CLOSE 

tt  i 

MOV 

<SP)f. R5 

MOV 

<SF\>+,  F'4 

MOV 

< SF‘  ;■  * . PS 

M 0 V 

<SP>+,  F:2 

MOV 

csp Ri 

MOV 

(SP >♦,  PO 

RTS 

PC 

BLKM 

5 

. RADSO.'DK  SHIHlODftT / 
. EVEN 
. CSECT 

. END 


FETCH,  . LOOKUP,  . F.'EADW 


, PUSH  PO 
. PUSH  ft 1 
: PUSH  R2 
, PUSH  RS 
; PUSH  ft 4 
; PUSH  ft 5 

■ SET  UP  Ri  WITH  HC'IPESS 
.SET  UP  ft?  ftS  BLOCK 
.•SET  UP  RS  ftS  DftTft  FTS 
.DEFINE  FILE 

FT  , PEfiD  FROM  DISK' 

■ CLOSE  CHANNEL  til 
, F OP  ft  5 

, F 0 F P 4 

■ p OF  R S 
, PDF  F 2 
. POP  PI 
, POP  RO 
, RETURN 

.DftTft  STORED  FiS  THIS  NOME 
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JO 

so 


50 


60 


70 

75 


S5 

AO 


J 20 

J 40 

1 60 

JSO 

?00 
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SUE  ROl'T  I NE  S10P1 

DIMENSION  K<J?n2>.  X<  500  > . V<500> 

COMMON  k , A , TIME..  IF  RE  01.  IFREQ2,  I AXIS,  ID1 
TVFE  SO 

r OPMAT  < IX.  ' THE  1ST  1 LOST  PTS.  ? THE  FREO  DEV.?') 
ACCEPT  50. JS1. JSC. IFOEV 
FORMAT (Sib) 

JS?»JS1*< JS2-JS1 > 2 2 
N = 0 

X v 0 > = 0 0 
T I M E * 0 
T I ME  1*0. 

00  AO  1=8. JS2 
IFOv(I  V EC!  CO  GOTO  ?0 
IFOu  I > NE  -100)  GOTO  60 

1 = 1*1 

I F<k<I>  GE  0>  GOTO  60 
T I M E 1 3 < 655  ■ ♦ ■ A&S(  \ < I > ) ' *80. 

1 = 1 + 1 

TIME-Kt  1 '•.■'1750  0 + TIME1  + TIME 
I F < I - J ? l > 85.70  ~5 
T I ME 2 = T I ME 
N = N + 1 

V<N  ■ =2  5000  C»2k \ I > 

L = N - 1 

X<N>*K< I >21250  +X(L > 

I F < I NE  JSS>  GOTO  85 
T I M E F * < X < N > + T I M E 2 > 2 8 0 . 

TIME  1=0 
CONTINUE 

T I ME *K < 2 > ♦ C C 3 ) 2 1 000  ♦ T I MEF 

LOOPED 

A = 0.  0 

E:«0  0 

GOTO  200 

TEST*2  5MFDEV 

GOTO  200 

TES  T » I Ft'EV 

GOTO  200 

TES T =0.  5+  I FDEV 

GOTO  200 

TE ST  =0  25* I FDEV 

NDAT A=0 

S U M X = 0 OEO 

SUM V*Q.  OEO 

SUMXX-0  OEO 

SUMXVaO  OEO 

DO  400  1*1, N 

I F (LOOP  EO  0>  GOTO  TOO 

F R E 0 = 0 25«<.  0 1 2 5 ♦ A ♦ X ( I > ♦ B > 

DEV  I *A8S «.  V v I >•  FRED  • 


! 


I 
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I F <DEVI GT.  TEST  > GOTO  400 
300  SUMX*SUMX*'XC  I > 

SUMV»SUMV  + V< I ) 

SUMXX**SUMXX*X •'  I >»  X<  I 
SUMXV  = SUMXV«-X\  I HV(  I > 

NDATA*N0ATA+1 
400  CONTINUE 

D = N D ft  T ft  * S U M X X - S U M X ♦ S U M X 
C = ND  A T A * S U M X V - 5 LI  M X »•  SUN  S' 

E = S U M V ♦ S U M X X - S U M X * S U M X V 

IFCNDATA.  L E 1.  OR.  0 E0.  O.  0)  GOTO  10 

A = 320  O + C/0 

0=4  O+E/O 

LOOP =100 F+l 

GOTO  < ISO. 140, 160. ISO. 50O>, LOOP 
500  SUMS=0  0 

DO  600  1=1, N 

FRE0  = 0.  25  ♦ ( 0 1 2 5 *■  A * X ; I > + B > 

D E V I = A B S ( V<  I > - F R'  E 0 > 

IF (DEVI  GT  TEST  > GOTO  600 
SUMS»SUMS  + DEV*I  **2 
600  CONTINUE 

5 U M S = S 0 R T < S U M S X N D A T A > 

TVPE  700,  NDA  T F . N,  A,  B,  SUMS,  TIME 
700  FORMAT  < IX.  ' 4» R m T I 0 # ' . 13,  ' /' , 13,  2X,  ' SLOP  : ' , FS».  3,  2X 

» ,' 6*' , F7.  3,  2X.  ' RES.  RMS  ' , F6.  3,  2X,  ' T I ME  : ' , F8.  3> 

I F < I C*  1 . NE  5 > GOTO  760 
GOTO  BOO 

760  I A = I N T (5  A , ' < IFRE02-IFPE01  > ) * ( I A X I S / 3 . 125)) 

1B  = INT<  <B-IFRE0i>*C32G0.  / ( I FREP2  - 1 FRE01 ) ) + 400.  > 
??0  CALL  DSLOPCK,  ID1,  JS1,  IA,  IB,  IAXIS) 

BOO  RETURN 

END 
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10 

30 

SO 


80 

?0 

75 


63 

90 

120 

140 

iso 

180 

200 
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SUBROUTINE  Si. OPS 

DIMENSION  K>.  17921.  Xt300'.  V<300> 

COMMON  k. A. TIME. IFREG1. IFREQ2, lANlS, ID1 
TVPE  30 

FORMAT < 1 X. 'THE  1ST  * LAST  PTS  ? THE  FREQ  DEV  ?' 1 
ACCEPT  50.  JS1.  JS2.  IFDEV 
FORMAT (313) 

JS3»JSl*v JS2-JS1 >/2 

N-0 

X<O>-0  0 

TIME-0 

TIMEl-P 

DO  90  1-8. JS2 

IF<K<  I >.  EP  0^  GOTO  90 

I F < k v I >.  NE  -1001  GOTO  80 

I-I  + l 

IFCKt  I > GE  O'1  GOTO  80 
T I ME  1 - < 853>*<.ABS(k<  l > > >*80 
I - 1 ♦! 

TIME-kCI V1230  ♦TIMEl+TIME 
IF<  WS1>  85, 7P, 75 
T I ME2-T I ME 
N"N*1 

V<N-> -25000.  0,‘k<  I > 

L-N-l 

X v N 1 ■ k < I ''.‘1250  ♦ XvL  > 

I F < I . NE  OSS'1  GOTO  85 
T I MEF  - 1 X<  N>  *T  I ME2  V80. 

TIME1-0 

CONTINUE 

TIME-KC'2  '♦k\  3>/l  000  ♦ T I M E F 
L OOP- 0 
A-0  0 
8-0  0 
GOTO  200 
TEST-2  5 ♦ I F P E V 
GOTO  200 
TEST- I FDFV 
GOTO  200 
TEST-0  5 * 1 F D E V 
GOTO  200 
TEST-0  25*If  DEV 
NDATA-0 
SUMX-0  OEO 
SUMV-0  OEO 
SUMXN-0  OEO 
SUMXV-0  PEP 
DO  400  I -1. N 
I F v LOOP  EP  P>  GOTO  300 
FREQ-0  23*  ( . P125*A*XlI>*f' 

0 E V I - A 6 S v V C I J-FREP) 
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IF<DEVI,GT.  TEST;-  GOTO  400 
:<oo  S U M X=SUMX+X<I) 

S U M V = 5 U M V + V C I :• 

sumxx«suhxx+ x < I >*>:c  I ) 

SUMXV=SUMXV+X< I)*V(I) 

NDATA=NDfiTfi*l 
40O  CONTINUE 

D = N 0 ft  T ft  + S U H X X - S U M X * 5 U H X 

C = N D ft  T ft  * £.  U M X V - S U M X * S U M V 

E = S U M V ♦SUMXX- S U H X + S U M X V 

I F < N D ft  T ft . LE.  1 OR.  0.  EG.  0.  d)  GOTO  10 

A = 3 £ 0 . 0 ♦ C / 0 

B = 4.  0+-E/D 

LOOP-LOOP+1 

GOTO  <120.  140.  160.-  100,  5 0 0 ;• ..  LOOP 
500  S U H 5 = 0.  0 

DO  bOO  1=1. N 

F R E G = 0 . 25 * < . 012 5 * ft ♦ X < I > ♦ 6 ) 

D E V I = ft  B S < V < I > - F R E 0 ) 

I F < D E V I . GT  TEST''  GOTO  60O 
SUMS=SUMS+DEVI**? 

600  CONTINUE 

S U M S = S 0 R T < S U M S / N D ft  T ft  ;• 

TYPE  TOO, NDftTh  N, ft, 6, SUMS, TIME 
700  FORMAT < 1 X.  ' #RAT 10#'  . l'<-  ' /'  , I i ■ 2 X • ' SL  OP  : ' , FS*.  3,  2X 

# , ' B = ' , FT.  Z , 2X,  ' RE S.  RMS  : ' , F6  3 . 2 X . ' TIME  ; ' . FS.  3 ) 
IFCIDl  NE  5'-  GOTO  76 O 
GOTO  r?o 

?60  I ft  = I N T < 5 . ♦<ft/< IFREG2-IFREGl))*<IflXIS/2.  125)) 

I B = I N T < < B - I F R E C 1 ) * <0200  /'  I FREG2- 1 FRFQ1  > )+400.  ) 
780  CftLL  DSLOPCK,  ID1. JSi,  Ift,  IB.  IftXIS) 

790  fi=-fi 

BOO  RETURN 

END 
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. TITLE 

V-AXIS  expansion 

G L 0 B L 

YEXPO. VSHOWi 

LPSVC*17 

'0410 

. HCALl 

V.1  •> 

. Vi  . REGOEF 

V «. 

REGDEF 

MACRO 

SHOW  A , B . C . 0 

, MACRO  CALL 

CMP 

Rl.HO 

, LOWER  LIMIT  >nO? 

BuT 

. + iii 

, VES. NOT  OISPLAV 

CMP 

ri.  no 

, NO. LOWER  LIMIT  CnO" 

BIT 

♦ ii 

, VES.  COMPARE  UPPER  LIMIT 

MOV 

Hi. FLAGVl 

. NO,  LOWER  LIMITMD 

BR 

. *40 

. GOTO  OISPLAV 

CMP 

Ri.  no 

. UPPER  LIMIT  OnO"' 

BGE 

♦ b 

. VES- KEEP  TRVING 

JMP 

FINISH 

, NO. F INISH ! 

CMP 

Ri.  no 

. UPPER  LIMIT 

BNE 

♦li 

, NO. GOTO  OISPLAV! 

MOV 

n.T.  FLAuVl 

, UPPER  LIMIT  «> 

BR 

♦ 10 

.GOTO  OISPLAV 

MOV 

ni, FLAGVl 

, no  IS  WITHIN  LIMITS 

MOV 

WH, VP0S1 

. SAVE  F I RS T C HARAC  TER 

MOV 

nB. VPOSi 

• SAVE  SECON'O  CHARACTER 

MOV 

nc . VFOSS 

.SAVE  THIRO  CHARACTER 

MOV 

no. VSCALE 

J SR 

PC . VSHOWi 

, SUBROUTINE  FOR  OISPLAV 

CLR 

FLAGVl 

EN[>M 

MOV 

RO.  -v.SP> 

, PUSH  RO 

MOV 

RI,  -tSF‘> 

. PUSH  Rl 

MOV 

Ri. -CSP> 

. PUSH  Ri 

MOV 

RS. -eSP' 

, PUSH  RS 

MOV 

R 4 , -v  SP'1 

. PUSH  R4 

MOV 

R5,  -(SP'' 

. PUSH  RS 

MOV 

n ioooo, i ps vc 

.ERASE  THE  SCOPE 

MOV 

i»i(RS'<,  FPE01 

; MOVE  LOWFR  LIMIT 

MOV 

P 4 v.  R 5 > . FREOi 

•MOVE  UPPER  LIMIT 

MOV 

0bCR5>. FLAGVi 

.VALUE  FOR  DIVISION  ON  SCOPE 

CLR 

FLAGVl 

. TEST  FLAG  FOR  nO 

MOV 

FRE01 . RI 

, SET  UP  Rl 

MOV 

FREOi. Ri 

, SET  UP  Ri 

SHOW 

SPACE • Nb.  NO-  £0 

, OISPLAV  BO  ANO  HONRI  LINE 

SHOW 

SPACE. NT- NO. 70 

.OISPLAV  70  ANO  HONRI  LINE 

SHOW 

SPACE. NS, NO. SO 

.OISPLAV  BO  ANO  HONRI  LINE 

SHOW 

SPACE.  NO.  Nu  ?0 

.OISPLAV  ?0  ANO  HONRI  LINE 

SHOW 

Nl, NO  NO. 100 

•OISPLAV  100  ANO  HONRI  LINE 

SHOW 

Nl.  Nl.  NO.  1 1 0. 

.OISPLAV  110  ANO  HONRI  LINE 

SHOW 

Nl . Ni. NO. liO 

.OISPLAV  liO  ANO  HONRI  LINE 

SHOW 

Nl  n:.no.  i:o 

•OISPLAV  ISO  ANO  HONRI  LINE 
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SHOW 

N 1 , N4  NO.  1 40 

SHOW 

Nl.  N5. NO,  150 

SHOW 

N 1 , N6,  NO.  160 

SHOW 

Nl,  NS,  NO,  ISO 

SHOW 

N2.  NO.  NO,  200 

FINISH 

MOV 

<SP>«-. R5 

MOV 

C S P ) ♦ , R 4 

MOV 

iSP'1*-.  PS 

MOV 

< SP  P2 

MOV 

V s P > ♦ • R 1 

MOV 

<rsp>*.  po 

RTS 

PC 

CSECT 

VDATA2 

VP0S1 

. ELI  W 

T 

VP0S2 ; 

. E'L  KW 

YP0S3 

EL  KM 

7 

VSCALE 

. WORD 

0 

FLuGYl 

WORD 

F LAGV2 : 

. WORD 

F RcOl  : 

WORD 

FREQ2 

. WORD 

. CSECT 

NLM'l 

NO 

. EYTE 

?e.  121,  in.  105-  7c 

N1 

BYTE 

O,  102,  1 77,  IOC,  0 

N2  : 

BYTE 

142. 121. 111. 105. 10 

NS  : 

. BYTE 

42, 101, 111. ill, 66 

N4 

. BYTE 

SO.  24.  22,  177,  20 

N5  . 

EYTE 

47,  105,  1 05,  1 05,  71 

No 

BYTE 

76, 111. 11 1, 111, 62 

N7 

BYTE 

101,  41,  21,11,7 

N3  . 

BYTE 

66,  111, 11 1. Ill  66 

N3  . 

BYTE 

46.  Ill,  11 1,  1 11 • 76 

SPAC  E 

EYTE 

0 , 0 • 0 . 0 , 0 

, [■  I SPLAY 
, DISPLAY 
, DISPLAY 
. DISPLAY 
. DISPLAY 
. POP  P5 
. POP  R4 
, POP  RS 
, POP  R2 
, POP  R 1 
, POP  Rd 
, PET  U R E 
, SECTION 


140  AND  HONRI. 
150  AND  HONRI 
160  AND  HONRI 
ISO  AND  HONRI 
200  AND  HONRI 
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LINE 

LINE 

LINE 

LINE 

LINE 


FOR  Y-AXIS  DISPLAY 


.AND  NUMERICAL  CHARACTERS 


,•  S E C T I 0 N (•  0 R N L ’ t.  V I i.  n L 


EVEN 

CSECT 

END 


YEXPD 
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C SECT 

T 1 TLE 

VS  HOMO 

GLOBl 

VS HO WO 

l PSVC  * 1 

ro-uc- 

l PSVCK« 

170400 

L PS  VC  S' = 

170422 

mcall 

V2  . . REGDEF 

V *. 

REGDEF 

•THIS  SUBROUTINE  DISPLAYS 

NUMERICAL  C HftRC  T E RS 

. ON  THE 

SC  OF  E ft  NO  EXPAND 

THE  VERTI  SCALE 

VSH0N2 

MOV 

F.  0 , - i.  S F > 

FUSH  RO 

MOV 

R 1 , - < S P s 

PUSH  R 1 

MOV 

R2.-v.SPs 

PUSH  R2 

MOV 

RS.-v.SP) 

F’  U S H R 

MOV 

R 4 . - < S P s 

PUSH  R 4 

MOV 

R5.-v.SF) 

PUSH  R5 

CL  R 

R4 

R 4 = 0 

CLP 

TEST 

TESTED 

CLR 

LFSVCX 

LPSVC X«0 

CMP 

• O.FLftGVl 

FLAG  S'  1=0? 

BNE 

ft  L P H ft  0 

NO. GOTO  OISPLftV 

J HP 

our 

VES. GO  OUT 

ALPHAO 

CMP 

PI. FlftGVl 

PLftGs'l=0*' 

BNE 

ALPHA 

NO. GOTO  ft L F H ft 

MOV 

#400  VL I NB 

VES- LOWER  LIMIT 

MOV 

VLINE.R5 

SAVE  R5 

BR 

GftMMft 

GOTO  GftMMft 

ALPHA 

CM^ 

#S.FLAGS'l 

FLAGY1=>G  •' 

BNE 

Ft  T ft 

NO.  GOTO  BET  ft 

MOV 

wrcoo  VLINE 

VES  UPPER  LIMIT 

MOV 

VLINE.F5 

SAVE  R 5 

Efi 

GftMMft 

GOTO  GftMMft 

PETR 

MOV 

VSCFUE  R1 

SET  UP  Pi  WITH  ADDRESS 

SUB 

F R c 0 1 • R 1 

Rl*R'l -LOWER  LIMIT 

C L F 

fi  5 

R 5 * 0 

DELTA 

add 

F L ft  b s'  2 . R 5 

R5  = SC  OF  E [■  I V I S I ON*  R 1 

DEC 

PI 

CMP 

PORI 

BLT 

OELTft 

hC  D 

#400  • R 5 

R5  = R5 ♦ 4 00  v BftSE  VALUE ' 

MCV 

P5 • VL INE 

SAVE  R5  FOR  HONRI  LINE 

GAMMA 

SUB 

#4  0.  R 5 

R5=-R5-;2 

MOV 

R 5 • S'  S T ft  R 

DISPLAV  1ST  CHARACTER 

MOV 

VF 0S1 ■ F 0 

BR 

MU 

GOTO  MU 

OMEGA 

MOV 

VP0S2.R0 

Of SF LAV  2ND  CHARACTER 

MOV 

VSTftR. R5 

DISPLAV  HONRI  LINE 

BR 

MU 

GOTO  MU 

THETA 

MOV 

VPOSS.RO 

DISPLAV  SRD  CHARACTER 

MOV 

S'  S T ft  R • R 5 

DISPLAV  HONRI  LINE 

Cl1  - NfW325.MAC,  Disolavlnu  SVrtlCrtl  Cpordln.it? 
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MU 

MOV 

# - 5,  R 1 

IOTA 

ADD 

#15,  R4 

MOV 

VS  TAR, R 5 

MOV 

#-?,  R 2 

MOVE: 

( R 0 > + , r: 

ETA  : 

ROLP 

R3 

PPL 

HU 

MOV 

#2002.  LF'S 

PHI  : 

TSTE 

L F'SVC 

PPL 

F HI 

MOV 

F 5,  L F'SVC  V 

MOV 

R 4 . LPSVCX 

NU  : 

ADD 

# 1 1 , R 5 

INC 

R2 

PNE 

ETA 

INC 

PI 

PNE 

IOTA 

INC 

TEST 

CMP 

#1. TEST 

PE  Q 

OMEGA 

CMP 

#2, TEST 

PEG 

THETA 

ADD 

#100.  R 4 

MOV 

YL  INE.  F'2 

SIGMA . 

MOV 

#2002.  IPS 

OK . 

TSTE 

L F'SVC 

PPL 

OK 

INC 

F.  4 

MOV 

F:  4 . LPSVCX 

MOV 

R2, LF  SVC  V 

CMP 

#7770, R4 

PGE 

SIGMA 

OUT  . 

MOV 

<SF  F:5 

MOV 

( S P '•  +■ ..  R 4 

MOV 

<■  sp r: 

MOV 

<SP  >+,  P2 

MOV 

< S P > + , P 1 

MOV 

< s p > + , r o 

RTS 

PC 

VST  AR . 

. WORD 

0 

YLINE 

. WORD 

0 

TEST. 

. WORD 

0 

. CSEC.T 

VDATA2 

YPOS1 : 

BLKW 

3 

VP0S2 : 

. BLKW 

3 

VP0S3 : 

. BLKW 

VSCALE . 

. WORD 

O 

FLAGVi . 

. WORD 

FLAGY2: 

. WORD 

F PE 01 

. WORD 

F RE 02 : 

. WORD 
. C SEC  T 

. END 

VSH0W2 

/ R 1 = - 5 FOR  5 COLUMNS 
/POSITION  FOR  NEXT  COLUMN 
/ INITIALIZE  V-POSI T I ON 
. P2  = -7  FOR  7 DOTS/COLUMN 
.SET  UP  R3  WITH  ADDRESS 
z ROTATE  R3 

zCRRRV  PIT  SET. DISPLAY  A DOT 
; INITIALIZE  THE  SCOPE 
. READ V ? 

. NO. WAIT 

j VES, MOVE  V-VALUE 
. VES.MOVE  X-VALUE 
.POSITION  FOR  NEXT  DOT 
.FINISH  7 DOTS? 

. NO. GOTO  ETA 
.FINISH  5 COLUMNS? 

. N 0 , G 0 T 0 IOTA 
.FINISH  A CHARACTER 

. TRY  END  CHARACTER 

■ TRY  3RD  CHARACTER 
;R4=P4+64  FOR  X-POSITION 

. INITIALIZE  THE  CSOPE 
. READY? 

■ NO.  WAIT 

. R4=R4+1  FOR  X-POSITION 
. FOR  HONRI.  LINE 
. THE  SAME  V-VALUE 
, 4 0 S 5 > = R 4 
j VES, GOTO  SIGMA 
, POP  R5 
. POP  R4 
, F 0 P R 3 
, POP  RZ 
, POP  R 1 
, POP  R0 
, RETURN 


; SECTION  FOP  V-VALUE  AND 
.NUMERICAL  CHARACTERS 
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. CSEC” 

. TITLE 

DISP.  THE  DATA 

. GLOBL 

XEXPD, XSH0W2 

L P S V C * 1 

70416 

L PS VCX* 

170420 

LPSVCV* 

170422 

•.  MCALL 

. . V2.  . , . REGDEF 

. . V2.  . 

. REGDEF 
. MACRO 

SHOW  A.B.C 

MOV 

# A.  XP0S1 

MOV 

#e, xpos2 

MOV 

#C, X P 0 S 3 

JSR 

PC, XSH0W2 

XEXPD : 

. END  M 
MOV 

PO, -<SP> 

MOV 

PI, -<SP> 

MOV 

R2,  - < S P > 

MOV 

R3,  -CSP> 

MOV 

R4,  -<SP> 

MOV 

R5,  -CSF  ) 

MOV 

2 < R 5 > , RO 

MOV 

0 4 < P 5 > , PI 

MOV 

06  ( F:5  > , R2 

MOV 

c- 1 o < p 5 , data:; 

MOV 

0 1 2 < R 5 1 » D A T A 4 

MOV 

914<R5>, XI 

MOV 

0 1 6 < R 5 > , VI 

MOV 

0 2 0 < R 5 > » LAST 

MOV 

Rl, DATA  5 

MOV 

R2. DATA2 

SUE 

Rl,  F:2 

MOV 

R2, DATA1 

CLR 

R 3 

DEVI  : 

SUB 

#3,  R2 

INC 

R3 

CMP 

#3,  R2 

BLE 

DEVI 

MOV 

R3 . 1 AST 

OMEGA 

ASL 

Rl 

ADD 

Rl,  RO 

DEC 

DAT  A5 

START 

CLR 

IPS vex 

CRT1  : 

CMP 

< RO  > , #-100 

BNE 

ALPHA 

ADD 

#4,  RO 

ADD 

#2, DATA? 

ADD 

#100.  , LF’SVCX 

MOV 

#3 POO  . R3 

JMP 

move: 

C20  - 

NEW335.MAC,  Displaying  Data 

OF  100K  HZ 


; MACRO  CALL  FOR  VERT  1 VALUES 


i PUSH  RO 
; PUSH  R1 
, PUSH  R2 
, PUSH  R3 
, PUSH  R4 
, PUSH  R5 

, SET  UP  RO  WITH  ADDRESS 
j 1 ST  DATA  POINT 
, LAST  DATA  POINT 
.LOWER  LIMIT 
.UPPER  LIMIT 
-DISPLAY  INDEX 
.DISPLAY  INDEX 
; INDEX  FOR  X-VALUE 
. R 1 = D A T A 5 
, R2-DATA2 
, R 2 * R 2 - R 1 

i SAVE  P2  AS  DATA  POINTS 

. r:=o 

, R3*R2/2 


. LAST*R2 
i R1*?*R1 

. RO* ADDRESS  OF  1ST  DATA  PT 
, F I R S T = F I R S T - 1 
, L P S V C X = 0 
; VALUE* -100? 

, NO. GOTO  Al PHA 
i SKIP  TWO  DATA  PTS 

> DISPLAY  IN  SPECIAL  FORM 

.GOTO  MOVES 


and  Horizontal  coordinate 


v 


\ 
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ALPHA 

MOV 

R0>*. R4 

. SAVE  OATA  IN  R4 

INC 

.[•AT  ft? 

. D A TA5 -D AT  A5 ♦ 1 

CMP 

DATA?.  DATA; 

, DATA5C-DATA2? 

BLE 

crt;- 

. VES,  GOTO  CRT!' 

JMP 

OUT 

, NO. GOTO  OUT 

CRT2 

MOV 

•;oo;. lpsvc 

. SET  UP  STATUS  OP  SCOPE 

READY 

TSTB 

LPSVC 

.SCOPE  RE  AOV? 

PPL 

RE  fit* S' 

. NO, WAIT 

CMP 

R4,  DATA: 

. PPEO  > LOWER  LIMIT' 

OLE 

GAMMA 

•VES. GOTO  GAMMA 

MOV 

#;50  . RT 

.NO. MOVE  250  TO  V-VALUE 

BR 

MOVE  1 

.GOTO  MOVE  1 

GAMMA 

CMP 

P 4 . DATA4 

. FREO  CUPPER  LIMIT"* 

BGE 

BELT  A 

. VES. GOTO  DELTA 

MOV 

ttssso  , f : 

.MOVE  SB50  TO  V-VALUE 

BR 

MOVE  1 

.GOTO  M0VE1 

f>  E L T fl 

CLR 

p; 

. p;«o 

CLR 

PS 

. PS-0 

MOV 

P4. PESEP1 

, SAVE  R 4 

CLR 

R 4 

, R 4 ■ 0 

MOV 

VI.  PI 

, SET  UF  R1 

THETA 

ROD 

pi.  r; 

. DOUBLE-PRECISION 

ft  DC 

PS 

, <r;>lr: >«i;bo*ri 

INC 

P4 

CMP 

ttlTSO  . P4 

BGE 

THETA 

MOV 

PESEF1 . P 4 

.SAVE  R4-N0.  OF  COUNTS 

CLR 

Fl 

, F 1=0 

ETA 

SUB 

R 4 , f ; 

, DOUBLE-PRECISION 

SBC 

PS 

, ri»lr;>cr:o/r4 

INC 

PI 

CMP 

•to,  FT 

BLT 

ETA 

CMP 

R 4 . f ; 

BLE 

F • ft 

SUB 

XI.  Pi 

, Ri-Pl-Xl 

MOV 

r i . r : 

, FT-Pl 

MOVE  1 

CLR 

FT 

, Rl-0 

MOV 

LAST. ft 

DEVS 

SUB 

F F 4 

; R1-P4  -’R; 

INC 

PI 

CMP 

FT,  P4 

Bl  E 

[•EVT 

ADD 

FT  , LPSVC  X 

.ADD  Ri  TO  X VALUE 

MOVES 

MOV 

F LPSVC  V 

. FT* V-VALUE  ON  SCOPE 

MOV 

F I , TEST 

MOV 

LPSVC  X,  XSCALE 

.SAVE  L PE VCX 

MOV 

[•A TPS.  FT 

CMP 

0 A T A 1 . HI  50 

• DATA  PTSM50? 

BGT 

CM  01 

. VES.  GOTO  OF  01 

CLP 

FT 
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: 

; 


r 


1 
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MOV 

Rl,  R4 

0K0 

sue 

#12, P4 

, R4/10-R2 

INC 

RI- 

CMP 

M12,  R4 

eiE 

OKO 

CMP 

MO.  R4 

, R 4 = 0 ? 

BED 

OKI 

, VES,  GOTO 

OKI 

JMP 

CRT1 

, NO. GOTO 

CRT1 

OKI  : 

sue 

#1?.  R2 

i Rr/io 

CMP 

Mir.  Rr 

BL  E 

0 k 1 

CMP 

MO.  Rr 

, Rr=o? 

BEQ 

01  01 

, VES. GOTO 

OKOl 

SHOW 

SPfiC  E, SPHCE, 

SPACE  .VERTI 

. LINE  ON 

SCOPE 

JMP 

CRT1 

, GOTO  CRT1 

OK  01 

CMP 

M 1 00  . Ri 

BNE 

Oku 

SHOW 

N 1 . NO,  NO 

. OISPLfiV 

100  fiNO 

VERTI. 

JMP 

CRT1 

; L I NE  ON 

SCOPE 

Oku 

CMP 

M2 00  . Rl 

BNE 

ok  r i 

SHOW 

Nr. NO, NO 

, OISPLfiV 

roo  fiNO 

VERTI. 

JMP 

f RT1 

.LINE  ON 

SCOPE 

0K?1 

CMP 

»:oo  , ri 

BNE 

0101 

SHOW 

n:, no. no 

, OISPLfiV 

: oo  RNO 

VERTI. 

JMP 

CPU 

■LINE  ON 

SCOPE 

OKU 

CMP 

M400  , Rl 

BNE 

01,41 

SHOW 

N4 , NO. NO 

. OISPLfiV 

400  fiNO 

VERTI 

JMP 

C RT1 

L INF  ON 

SCOPE 

0 K *1 1 ; 

CMP 

M500.  . Rl 

BNE 

OK  51 

SHOW 

N5,  NO,  NO 

, OISPLfiV 

500  fiNO 

VERTI 

JMP 

CPT1 

.LINE  ON 

SCOPE 

OK  51 

CMP 

M600.  . Rl 

BNE 

OK  El 

SHOW 

NS.  NO.  NO 

i OISPLfiV 

BOO  fiNO 

VERTI 

JMP 

CRT  1 

; LINE  ON 

SCOPE 

OK«Jl 

CMP 

Mr 00  , Rl 

BNE 

ch  ri 

SHOW 

NT.  NO,  NO 

• OISPl fiv 

TOO  fiNO 

VERTI 

JMP 

C P T 1 

•LINE  ON 

SCOPE 

Ok?l 

CMP 

MBOO  . Rl 

BNE 

01  SI 

SHOW 

NS.  NO,  NO 

, OISPLfiV 

BOO  RNO 

VERTI 

JMP 

CRT1 

LINE  ON 

SCOPE 

OKS  1 

CMP 

MBOO  , Rl 

BNE 

0KB1 

SHOW 

NB,  NO. NO 

, OISPLfiV 

BOO  fiNO 

VERTI 

JMP 

C-PT1 

-LINE  ON 

SCOPE 

C 20  ~ NEW335.MAC  (Conti nurd) 
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0K31 : 

CMP 

#1000  , Ri 

ONE 

Oh  101 

SHOW 

Ml,  NO. NO 

i DISPLAV 

100  AND 

VFRTl 

JMP 

CRT1 

.LINE  ON 

SCOPE 

OK  1 01  ; 

CMP 

#1*100  , Rl 

ONE 

Ok  ill 

SHOW 

Nl.  N4,  NO 

• DISPLAV 

140  find 

VERTI 

OKI  11 

JMP 

CRM 

LINE  ON 

SCOPE 

OUT 

MOV 

LOST,  »20<R5> 

MOV 

v SP>  R5 

. PDF  R5 

MOV 

<SP,'f.  R4 

. F OF  R4 

MOV 

kSP  *f . R'S 

. POP  R2 

MOV 

<SP,<  + . R2 

. POP  Ri 

MOV 

t SP)  + . Rl 

. POP  Rl 

MOV 

CSP>  + , RO 

. POP  RO 

RTS 

FT 

. RETURN 

DATA!  ; 

WORD 

0 

DATA2. 

. WORD 

DATAS 

WORD 

DA  TR4 

. WORD 

DATA5 

. WORD 

0 

XI  . 

WORD 

VI 

. WORD 

RE  SERI 

WORD 

RESER2: 

. WORD 

LOST 

WORD 

CSECT 

NDATR2 

. SECTION 

FOR  X- VALUE 

XP0S1  : 

. PLkW 

XF0S2 

. PL  kW 

XP0S3  . 

. PLkW 

j*. 

XSCALE  . 

. NORD 

0 

TEST 

. WORD 
CSECT 

NO  Ml 

. SECTION 

FOR  NUMERICAL 

NO. 

. IVTI 

76,  121.  111.  103. 

7 e>  i 

CHARACTERS 

Ni 

. PVTE 

o.  102.  i 77.  ioo.  o 

N2  . 

. PVTE 

142.  121.  111.  105 

102 

N3 

. PVT F 

42. 101  111111. 

bb 

N4 

. PVTF 

SO.  24.  22.  172.  20 

N3 ; 

. PVTE 

47.  105.  105.  105. 

71 

Nb : 

. PVTE 

7b, 111. 111.  111. 

f>2 

NT : 

. PVTE 

101,  41, 21.  11,  7 

N3; 

PVTE 

t'b.  Ill,  111.  Ill, 

t*b 

N? 

PVTF 

4b, 111.  ill.  Ill, 

7b 

SPRCE  . 

. PVTE 
. EVEN 
CSFCT 

0.  0.  0 0.  0 

. END 

XEXPD 

P20  - NEW335.MAC  (Continuous 


. CSCCT 
TITLF 

MS HO M2 

. GLOBL 

XSH0IJ2 

L PSVC«i 704  IS 

i P5VCX- 

170410 

LPSVCV* 

170422 

. MCALL 

. . V2  . 

REGDEF 

V2  , . REGDEF 

i 

, THIS  S 

UBPOUTINE  IS  TO 

DISPLAV  NUMBER I CAL  CHARCTERS 

, ON  THE 

SCOPE  AND  EXPAND  THE  HONRI.  SCALE 

XSH0U2 : 

MOV 

RO. -<SP) 

, PUSH  RO 

MOV 

PI,  -<5F> 

; PUSH  Rl 

MOV 

P 2 , - <SP> 

» RUSH  R2 

MOV 

RZ. -<SP> 

, PUSH  RZ 

MOV 

R4. -CSR> 

i PUSH  R4 

MOV 

R 5 . -<SP > 

, PUSH  R5 

Cl  R 

R4 

, R 4 = 0 

CLR 

L PS VCX 

, L P S V C X = 0 

Cl  R 

LPSVCV 

, L P S V C V = 0 

MOV 

X SCALE ■ XL  INF 

.SAVE  XSCALE 

SUB 

#2t . , X SCALE 

; XSC  RLE=XSCALE-Z6 

MOV 

XP0S1 , RO 

i FIRST  CHARACTER 

BR 

GAMMA 

, GOTO  GAMMA 

ALPHA  . 

MOV 

XP0S2-  RO 

, 2ND  CHARACTER 

BR 

GAMMA 

; GOTO  GAMMA 

BETA: 

MOV 

XPOSZ, RO 

, CRD  CHARACTER 

GAMMA . 

MOV 

# - 5 • R 1 

, R 1 = - 5 FOR  5 COLUMNS, ’CHARAC. 

DELTA. 

ADD 

#14,  XSCALE 

•POSITION  FOR  NEXT  COLUMN 

CLR 

R5 

; R 5 = 0 

MOV 

#-7,  P2 

• R2=-7  FOR  7 ROWS/CHARAC. 

MOVE 

(.  R 0 > ♦ , R ’■ 

.SAVE  DATA  IN  RZ 

ZETA  . 

ROLE 

RZ 

; ROTATE  RZ 

BPL 

RAP  PA 

, CARRV  BIT  SET, DISPLAV  A DOT 

MOV 

#2002, L PSVC 

> SE  T UP  SCOPE 

PHI  ; 

T STB 

LFSVC 

, SCOPE  READV? 

BPL 

PHI 

; NO,  WAIT 

MOV 

R5. LPSVCV 

iVES, PUT  A DOT  ON  SCOPE 

MOV 

XSCALE, L PS VCX 

» APPA . 

ADD 

#11. R5 

V-F  OS  I T I ON  FOR  NEXT  DOT 

INC 

R2 

; R2=R2*  1 

BNE 

ZETA 

•FINISH  A COLUMN? 

INC 

Ri 

, R 1 = R 1 ♦ 1 

BNE 

DELTA 

.FINISH  A CHARACTER? 

INC 

R 4 

, F-4.P4  + 1 

CMP 

# 1 . R 4 

, P 4 3 1 ? 

BED 

ALPHA 

• VES.TRV  2ND  CHARACTER 

CMP 

#2,  R4 

# P 4 3 2 ^ 

BED 

BETA 

.VES.TRV  ZRD  CHARACTER. 
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. CSECT 
. TITLE 

DSLOP 

; DISPLAY  ft  FITTED  STRftIGHT 

. GLOBL 

DSL  OP 

.LINE  ON  SCOPE 

L P S V C = 1 

70416 

LPSVCX* 

170420 

L P S V C V 3 

170422 

. MCftLL 

. . V?.  . 

. REGDEF 

. . V2  . , . REGDEF 

DSL  OP . 

MOV 

RO, -<SP> 

, PUSH  RO 

MOV 

Rl. -<SP> 

. PUSH  Rl 

MOV 

R2, -<SP> 

, PUSH  R2 

MOV 

R3» -<SP> 

, PUSH  R‘ 

MOV 

R4 , -<SP> 

, PUSH  R4 

MOV 

R5» -<SP> 

, PUSH  R5 

MOV 

2<P5>, RO 

•SET  UP  RO  WITH  ADDRESS 

MOV 

C*4(R5>.  FIRST 

, FIRST =5 

MOV 

(?  6 < . R 5 > » R 4 

.SET  UP  1ST  DATA  POINT 

MOV 

C*  1 0 <.  R 5 1 > I A 

, I A = S L 0 P E 

MOV 

012(R5>,  IB 

, I B- INTERCEPT 

MOV 

C-l  4 < R5>,  LAST 

, LAST=DISPLAV  INDEX 

CLR 

R3 

. p;  .-o 

SUB 

FIRST,  R4 

, R 4 = R 4 - 5 

MOV 

FIRST, R2 

. P2=F I PST 

asl 

R2 

, F:  2 3 R 2 * 2 

ROD 

R2,  RO 

, ADDRESS  OF  1ST  DATA  PT 

MOV 

LOST, R2 

• SAVE  R2 

OKO0  : 

MOV 

i R O ) ♦ , R 1 

. MOVE  DATA  TO  Rl 

OK  01  . 

SUB 

R2.  Rl 

; R 1 . ‘ R 2 e R 3 

INC 

r: 

CMP 

R2,  Rl 

BLE 

OK  01 

DEC 

R4 

. R4*R4-1 

CMP 

#0.  R 4 

i 0 3 R 4 ? 

BLE 

Ok  00 

, VES. KEEP  DOING 

MOV 

r:, vex 

, X-VftLUE  OF  1ST  DATA  PT 

CLR 

LPSVCX 

, L PVCX«0 

CLR 

LPSVCV 

. L P V C V = 0 

Cl  R 

RO 

. R 0 = 0 

MOV 

ve  x,  R4 

, SAVE  R4 

MOV 

I ft.  R* 

, R 3 = I A 

MOV 

IB,  R2 

. R 2 * I B 

OK 3 5 : 

MOV 

#2002. LFSVC 

SET  UP  THE  SCOPE 

Ok  4: 

TSTB 

L PSVC 

. SCOPE  REA D V ? 

BPL 

Ok  4 

. NO,  WAIT 

CMP 

#0.  PO 

. VES,  R0  = 0"' 

BNE 

BET  ft 

, NO. GOTO  BETA 

sue 

R3.  R2 

, VES. R2-R2-R3 

sue 

#50  . R4 

, R4=R4-50 

CMP 

#0.  R2 

, 0>«R2? 

l 

\ 
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BGE 

GAMMA 

. V E S , GOTO  GAMMA 

CMP 

#0,  R 4 

. 0 < R 4 ? 

BLT 

OMEGA 

; V E S . GOTO  OMAGFt 

GAMMA . 

MOV 

VCX,  R 4 

; SAVE  X-POS I T I ON 

MOV 

I A,  PS 

.SAVE  SLOPE 

MOV 

IB.  R2 

.SAVE  INTERCEPT 

MOV 

# 1 . R 0 

. RC'  = 1 

BR 

Ok  4 

. GOTO  OK 4 

BETA 

a or* 

RSi  R2 

, R2«R2*R3 

ADO 

#50.  . R4 

> R4»R4+50 

CMP 

#rrre, R4 

. 40950  R 4? 

BL  E 

OUT 

. VES.GOTO  OUT 

CMP 

#r?rc, rc 

. 4095OP2? 

BLE 

OUT 

.VES.GOTO  OUT 

OMEGA . 

MOV 

R4. L PS VCX 

jDISPLAV  A DOT 

MOV 

R2.  LPSVCV 

BR 

OK  4 

j GOTO  Ok 4 

OUT  : 

MOV 

<SP>*.  R5 

, POP  R5 

MOV 

<SP>+.  R4 

. POP  R4 

MOV 

<SP>+. R2 

, POP  R3 

MOV 

. S P > ♦ » R'2 

, POP  R2 

MOV 

<SP>*. Ri 

. POP  Rl 

MOV 

< SF> ♦. RO 

. PDF  RO 

RTS 

PC 

. RETURN 

VCX  : 

. WORD 

0 

I A 

WORD 

0 

IB  ; 

. WORD 

0 

r I R5  T ; 

WORD 

LAST  : 

. WORD 

ENT 
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DISP 


DEVI 


OMEGA 


START 

CRT1 


ALPHA 


C5ECT- 


. TITLE 

DISP  THE  DATA 

OF  100k  H2 

. GL-OE'L 

[■ISP.  XSHOW 

LPSVC =1 

704  IE 

lfsvcx* 

1704  20 

L P S V C V = 

170422 

. 11  CALL 

r* 

. V 2 , . REGDEF 

v c . 

. REGDEF 

. MACRO 

SHOW  A.E'.C 

.MACRO  CALL  FOR 

MOV 

ftA.  XP0S1 

, DISPIAVING  X-AXIS 

MOV 

#B.  XP 0S2 

MOV 

ft  C . XPOSS 

JSR 

PC • XSHOW 

E NDM 

MOV 

PO. -<SP> 

, PUSH  PO 

MOV 

PI.  -•  SP> 

. PUSH  Rl 

MOV 

P2.  -i  SP  ■ 

, PUSH  R2 

MOV 

PC. -» ST  ' 

, PUSH  R*. 

MOV 

P4 . SP> 

, PUSH  R 4 

MOV 

F5,  -c  SP  . 

• PUSH  R5 

MOV 

2 ■ R 5 •'  PO 

.SET  UP  RO  WITH  ADDRESS 

MOV 

0 4 < R 5 > . R 1 

.SET  UP  Rl 

MOV 

C°  E ( F'  5 '•  ■ P 2 

.SET  UP  R2 

MOV 

i?10'  R5  ;■ . s-eSERl 

• SET  UP  P E SERI 

MOV 

PI  F IPST 

• SAVE  Ri 

MOV 

P 2. DATA  2 

, SAVE  R2 

SUE 

PI.  P2 

• F2=R2-R1=DATA  POINTS 

MOV 

R2.  DAT  Hi 

, SAVE  NO.  OF  DATA  PTS 

CLP 

P5 

. PS  = 0 

S 'J  E 

#4,  R2 

• R5  = R2V4 

INC 

R5 

CMP 

(t  4 P2 

ELE 

DEVI 

MOV 

F 5,  LAST 

. SAVE  R5 

flSL 

Ri 

, F 1 = 2 ♦ R 1 

ADD 

PI  • FO 

• F'0  = ADDRESS  OF  FIRST  DAT 

MOV 

FIRST. PI 

;R1*N0  OF  FIRST  DATA 

DEC 

PI 

CLP 

LPSvc 

. LPSVC X = 0 

CMP 

CFO’',  ft  -l  uO 

. VALUE  = -100'’ 

BNE 

ALPHA 

•NO. GOTO  ALPHA 

ADD 

ft  4 . PO 

.SKIP  TWO  DATA  PTS 

ADD 

ft  2 . P 1 

ADD 

t!  100  , LPSVC  A 

L PS  VCX* L PS VCX «■  iOO 

MOV 

ft  I $00  . f 5 

. DISPLAV  IN  SPECIAL  FORM 

JNF 

move: 

MOV 

<.  R 0 > » . F 4 

, SAVE  DATA  IN  R4 

INC 

PI 

. Rl =F1 *1 

CMF 

Pi.  DAT A2 

. pi :data2- 

i 


4 
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GflMMfl2 


DEV2  . 


MOVE 1 ; 
DEV 3 : 

MOVE 3. 


APPENDIX  C 

BLT 

CRT2 

VES, GOTO  CRT2 

JMP 

OUT 

GOTO  OUT 

CRT2  : 

MOV 

#2002, LPS VC 

SET  UP  THE  SCOPE 

REftOV; 

.TSTB 

LPS  VC 

SCOPE  REfiDV? 

BPL 

R E ft  D V 

NO,  WftIT 

CMP 

#480.  , R4 

480CN0.  OF  COUNTS? 

BLT 

G ft  M M ft  0 

VES, GOTO  G ft M M ft 0 

MOV 

#2500.  , R5 

DISPLfiV  IN  SPECIAL 

BR 

M0VE1 

GOTO  M0VE1 

GflMMftG : 

CMP 

#0,  R4 

0CR4? 

BLT 

G ft  M M ft  1 

VES, GOTO  GftMMftl 

MOV 

#200.  , R5 

R 5 = 200 

ADD 

#50.  , LPSVCX 

JMP 

MOVES 

GOTO  MOVES 

GflMMRl  : 

CMP 

#2000.  , ft 4 

BUT 

GfiMMfi2 

MOV 

#330.  , R5 

#7,  lps vc.': 

MOVES 
#65000, S2 
#30,  r:< 

R5 

R4,  Rt 

R3 

R5 

#0,  R3 
DEV  2 
R4,  F:2 
DEV2 
R3 

LAST# R2 
R 2 , R 4 
RS 

P2,  R4 
DEVS 

RS# LPS VCX 
P5, LPSVCV 
R5» TEST 
LPS  VCX,  X SC  RLE 
0 ft  T ft  1 , #150. 

OK  01 
RS 

Rl,  R4 
#12,  R4 
RS 

#12,  R4 
OKO 
#0,  R4 
OKI 
CRT1 
#12,  RS 


, DOUBLE-PREC I S I ON 


, R5=l, 600. 00O/R4 


; RS  = 0 

; SET  UP  R2 


, DISPLAY  ft  DOT  ON  SCOPE 


i D ft  T ft  P T S > 150? 
<VES,  Goto  0 K 0 1 


, R 4 / 1 0 = ft  t 


, R 4 = 0 ? 

# VES, GOTO  OKI 
, NO, GOTO  CRT1 
# RS/10 
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•; 


CMP 

#1 r: 

OLE 

OKI 

CMP 

#0,  PS 

, p:=o? 

BED 

OK  01 

. VES, GOTO  0 K 0 1 

SHOW 

S F A C E , SPACE. 

SPACE  ; VERTI.  LINE  l 

ON  SCOPE 

JMP 

CPT1 

; GOTO  CPT1 

OK  01  : 

CMP 

#100.  , PI 

ONE 

0111 

SHOW 

Nl,  NO  NO 

/MSPLftV  100  ft  NO 

VERTI. 

J M P 

CPT1 

■LINE  ON  SCOPE 

OKU 

CMP 

#200  - PI 

ONE 

ok  ?i 

SHOW 

N2.  NO.  NO 

• MSPLftV  200  AND 

VERTI. 

JMP 

CPU 

•LINE  ON  SCOPE 

0 K 2 1 : 

CMP 

#:oo  / pi 

BNE 

OK*  1 

SHOW 

NT,  NO/  NO 

/MSPLftV  SCO  AN D 

VERTI 

JMF 

CPU 

■LINE  ON  SCOPE 

ck:i  . 

CMP 

#400  , PI 

BNE 

0 K 4 j 

SHOW 

N 4 / NO. NO 

• MSPLftV  400  HNC' 

VERTI . 

JMF 

C F'Tl 

■LINE  ON  SCOPE 

OK  41 

CMP 

#500  • PI 

BNE 

CM  51 

SHOW 

N 5 . N 0 . N 0 

/MSPLftV  500  ft N D 

VERTI 

JMP 

C F T 1 

■LINE  ON  SCOPE 

OK  51  . 

CMP 

#800.  , PI 

BNE 

01  81 

SHOW 

NS,  NO.  NO 

, DISFLAV  BOO  AND 

VERTI. 

JMF 

CPU 

; l INE  ON  SCOPE 

0 K 6 1 

CMP 

#700.  / PI 

BNE 

CM  71 

SHOW 

N7,  NO,  NO 

■ MSFLAV  700  AND 

VERTI. 

JMP 

CRT1 

•LINE  ON  SCOPE 

01  71 

CMP 

#800.  , PI 

BNE 

OK  81 

SHOW 

NC,  NO.  NO 

, MSPLAV  BOO  AND 

VERTI. 

JMP 

CPT1 

•LINE  ON  SCOPE 

OK  81  ; 

CMP 

#800.  , PI 

BNE 

OK  81 

SHOW 

N8,  NO,  NO 

■ MSFLAV  800  AND 

VERTI. 

JMP 

CF  :T1 

•LINE  ON  SCOPE 

Ot  81 

CMP 

#1000  . F:l 

BNE 

CM.  101 

SHOW 

til.  NO,  NO 

.MSPLAV  100  AND 

VERTI 

OK  101 

JMP 

CPT1 

•LINE  ON  SCOPE 

OUT 

MOV 

LOST,  F.ESFP1 

•SAVE  PESEP1 

MOV 

( S F’  > ♦ . S'  5 

■ POP  P5 

MOV 

, S F' > ♦ , R 4 

. FOP  P 4 

MOV 

•:  S P , F'3 

, F 0 P P 3 

MOV' 

S F "■  + , F :• 

■ F OF  P2 

C23  - NEW333.MAC  (Continued) 

V 


i 


I 


1 


APPENDIX  C 


MOV 

S P > ♦ . Pi 

MOV 

( S F 1 ♦ . P'  0 

RTS 

PC 

PESER1  . 

. MQRE> 

DATA1 

. WORD 

0 

DATA2 . 

WORD 

. CSECT 

SCOFF 

FIRST  : 

WORD 

LAST 

. WORT* 

CSECT 

X DATA 

XP0S1 

. BLkW 

j. 

XPQS2 : 

BLKW 

“■ 

XP0S3 

EL  KM 

3 

XSCALE 

MORt> 

0 

TEST  . 

WORD 

CSECT 

NUMBER 

NO 

B V T E 

re,  i2i.  iii. 

N 1 

BVTE 

o lor.  i rr.  i 

N2 

E-VTE 

142.  121  • ill 

N3 

BVTE 

42. 101 , 111 , 

N4 

BVTE 

-■  0 . i 4 ■ t c.  . 1 . 

N5 

BVTE 

42.  105.  105- 

NS 

BVTE 

re  lii.iii, 

NT 

BVTE 

101.41.21.1 

NS 

BVTE 

66, 11 . ■ 111. 

M? 

BVTE 

4 6 . 11..  111. 

SPACE 

. BVTE 
EVEN 
. CSECT 

0.0  0 o.o 

END 

CISP 

, POP  Pi 
, POP  RO 
. RETURN 


» SECTION  FOP  DATA 


• SECTION  FOP  X- VALUE 


•SECTION  FOP  NUMERICAL 

10  5..  -£  .CHARACTERS 

00.  0 

. 105  • J 0c. 

ill,  E r'. 

r . i 

105-  1 

11  i ■ c : 
i , r 

in.  - f 

■111.  'S 
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